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The determination of the energy spectra of small spin systems as for instance given by magnetic 
molecules is a demanding numerical problem. In this work we review numerical approaches to 
diagonalize the Heisenberg Hamiltonian that employ symmetries; in particular we focus on the spin- 
rotational symmetry SU{2) in combination with point-group symmetries. With these methods one 
is able to block-diagonalize the Hamiltonian and thus to treat spin systems of unprecedented size. 
In addition it provides a spectroscopic labeling by irreducible representations that is helpful when 
interpreting transitions induced by Electron Paramagnetic Resonance (EPR), Nuclear Magnetic 
Resonance (NMR) or Inelastic Neutron Scattering (INS). It is our aim to provide the reader with 
detailed knowledge on how to set up such a diagonalization scheme. 
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I. INTRODUCTION 

Magnetism is a research field that is almost as old as 
human writing. It took several thousand years until its 
nature, which is quantum, could be determined. In 1928 
Werner Heisenberg published his work on the theory of 
ferromagnetism Zur Theorie des Ferromagnetismus^ in 
which he introduced what is today called the Heisen- 
berg model. ^ That this spin-only model is successfully 
applicable to magnetism rests on the property of many 
iron group elements to possess a quenched angular mo- 
mentum in chemical compounds.^ Therefore, for many 
magnetic substances the Heisenberg Hamiltonian pro- 
vides the dominant term whereas effects connected to 
spin-orbit interaction are treated perturbatively in these 
systems. For theoretical work on non-Heisenberg systems 
see e.g. Refs. SHZl 

The Heisenberg Hamiltonian 

^Heisenberg = - ^ Jij^{^) ' (1) 

models the magnetic system by a sum of pairwise interac- 
tions between spins. The interaction strength (exchange 
parameter) between spins at sites i and j is given by a 
number Jij with Jij < referring to an antiferromag- 
netic and Jij > to a ferromagnetic coupling. The spins 
are described by vector operators. 

In order to understand magnetic observables such 
as magnetization, susceptibility, heat capacity or EPR, 
NMR and INS spectra the knowledge of the full energy 
spectrum of the investigated small magnetic system as 
for instance a magnetic molecule is often indispensable. 
Although the Heisenberg Hamiltonian, Eq. ([T]), appears 
to be not too complicated, analytical solutions are known 
only for very small numbers of spinJ^KSI or for instance 
for the spin- 1/2 chain via the Bethe ansatz.^^ The at- 
tempt to diagonalize the Hamilton matrix numerically is 
very often severely restricted due to the huge dimension 
of the underlying Hilbert space. For a magnetic system 



of N spins of spin quantum number s the dimension is 
(25 + 1)^ which grows exponentially with TV. 

Group theoretical methods can help to ease this numer- 
ical problem. A further benefit is given by the characteri- 
zation of the obtained energy levels by quantum numbers 
and the classification according to irreducible representa- 
tions. This review intends to provide an overview of the 
latest developments in efficient numerical diagonalization 
techniques of the Heisenberg model using symmetries. 
In particular we focus on the spin-rotational symmetry 
SU (2) in combination with point-group symmetries. 

The full rotational symmetry of angular momenta has 
been employed for quite a while. In quantum chem- 
istry the method of irreducible tensor operators was 
adapted to few spin systems along with the upcoming 
field of molecular magnet ism. ^^^H^ Nowadays the com- 
puter program MAGPACK^ that completely diagonalizes 
the Heisenberg Hamiltonian using SU(2) symmetry, is 
freely available. Also for the approximate determina- 
tion of energy eigenvalues by means of Density Matrix 
Renormalization Group (DMRG) methods, that can 
for instance treat chains of a few hundreds of spins with 
high accuracy, spin-rotational symmetry was employed.!^ 
In other fields such as nuclear physics this method was 
also adapted to model finite Fermi systems such as nu- 
clei employing SU{2) symmetry.^ Early applications are 
also known for Hubbard models, where one can actually 
exploit two SU{2) symmetries .231^ 

Besides spin-rotational symmetry many magnetic 
molecules or spin lattices possess spatial symmetries that 
can be expressed as point-group symmetries. Neverthe- 
less, a combination of SU{2) with point-group symme- 
tries is not very common. The reason, as will become 
more apparent later, is that a rearrangement of spins due 
to point-group operations easily leads to complicated ba- 
sis transformations between different coupling schemes. 
A possible compromise is to use only part of the spin- 
rotational symmetry (namely rotations about the z-axis) 
together with point-group symmetrie^^®! to expand 
all basis states in terms of simpler product states.'^^^*^ 
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During the past years a few attempts have been un- 
dertaken to combine the fuh spin-rotational symmetry 
with point-group symmetries. Ohver Waldmann com- 
bines the fuh spin-rotational symmetry with those point- 
group symmetries that are compatible with the spin cou- 
pling scheme, i.e. avoid complicated basis transforma- 
tions between different coupling schemes. Along the 
same lines, especially low-symmetry groups such as D2 
are often applicable since the coupling scheme can be 
organized accordingly, compare Ref. |35| for an early in- 
vestigation. Sinitsyn, Bostrem, and Ovchinnikov follow 
a similar route for the square lattice antiferromagnet by 
employing point-group symmetry. 

Very recently a general scheme was developed that al- 
lows to combine spin-rotational symmetry with general 
point-group symmetries.^^ The key problem, that the 
application of point-group operations leads to states be- 
longing to a basis characterized by a different coupling 
scheme whose representation in the original basis is not 
(easily) known, can be solved by means of graph the- 
oretical methods that have been developed in another 
context .f^^J^ We discuss in detail how this method can 
be implemented and present results for numerical exact 
diagonaliza tions of H eisenberg spin systems of unprece- 
dented size) ^^ * ^^ * ^^ * ^^^ Our aim is to provide the reader 
with sufficient material to be able to employ these pow- 
erful group theoretical methods. They can for instance 
also be applied to calculate higher order Wigner-n J sym- 
bols that appear when the double exchange is modeled 
in mixed-valent spin systems. ^^ ^^ 

The article is organized as follows. Section [ll] intro- 
duces the basic concepts, i.e. the Hamiltonian and its 
properties, the irreducible tensor operators, point group 
operations, and the construction of basis states for irre- 



ducible representations. Section [nl] demonstrates with 
the help of three examples that the Hamiltonian of spin 
systems of unprecedented size can be diagonalized com- 
pletely. The outlook in Section |IV| shortly summarizes 
and shows perspectives. The main part of this article 
is contained in an extended appendix that explains all 
technical details to set up the discussed diagonalization 
scheme. 



II. CONCEPTUAL IDEAS 

A. Spin Hamiltonian of magnetic molecules 

The research field of molecular magnetism deals with 
the investigation of the magnetic properties of chemical 
compounds composed of a number of ions that reaches 
from only a few up to dozens of it. For the magnetic 
modeling of the molecule only those ions are taken into 
account which possess unpaired electrons and thus a non- 
vanishing magnetic moment. Since the molecules, which 
are prepared in the form of a crystal or powder sample, 
are often quite well separated from each other by their 
ligands, inter-molecular interactions can be neglected in 



most cases. Additionally, the electrons can very often be 
treated as localized so that both features lead to a sim- 
plified sketch of the chemical compound, namely a spin 
system. The interactions between different spins of the 
system then depend of course on the chemical surround- 
ing and stem from direct exchange or super-exchangd^ 
via chemical bridges. Figure [l] shows such a simplifica- 
tion for a Crg compound which can easily be modeled by 
a ring-like system of interacting spins .l^^^ 





FIG. 1: Simplification of the chemical structure of a Org com- 
pound (l.h.s.) to the corresponding spin system for the same 
molecule (r.h.s.). Dots represent spin sites, lines exchange 
interactions. 



A general Hamiltonian that models magnetic molecules 
could be written as 



general 



H. 



exchange 



H 



Zeeman 



(2) 



where, to be more specific, in a system of N spins, i.e. 
z, j = 1, . . . , A^, the two terms assume the form 

-^exchange ~ ^ ^ S(^) ' ^'^j ' S^*^) ' 

^Zeeman = l-^B^ S (z) - gi ■ B . (4) 

i 

^exchange dcscribcs in a compact way the (isotropic and 
anisotropic) exchange interaction between two single-spin 
vector operators s{i) and s(j) as well as the single- ion 

anisotropy.^^ The quantity Jij is a second rank Carte- 
sian tensor containing the corresponding parameters. 
^Zeeman couplcs the spius to an external magnetic field 

B. In general, the coupling to an external field can be 
anisotropic and is thus described by local tensors g^. 

It turns out that for many magnetic molecules the 
isotropic Heisenberg Hamilton operator provides a very 
good model. In addition, we assume that for the highly 
symmetric spin systems to be treated in this article the 
^-tensors are scalars, and the same for all ions. Then the 
resulting Hamiltonian that models the system simplifies 
to 



H = H 



Heisenberg 



Zeeman 



(5) 
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with 

^Heisenberg = - ^ Jij s{i) - s{j) (6) 
Hzeeman = 9I^B^ ' B . (7) 

S = s{i) is the total spin. As already mentioned, 
Jij < refers to an antiferromagnetic and Jij > to a 
ferromagnetic coupling. 

The Heisenberg Hamiltonian is completely isotropic in 
spin space {SU{2) symmetry), i.e. the commutators of 
the square of the total spin S and its z-component 

with ^Heisenberg Vauish 



H 



Since 



Heisenberg? 







^Heisenberg, §^ — . (8) 



= the total magnetic quantum number 

M and the quantum number of the total spin S serve as 
good quantum numbers and a simultaneous eigenbasis of 

and ^Heisenberg CaU be fouud. 

A well adapted basis is then given by states of the 
form \a S M) which can be constructed according to a 
vector-coupling scheme (see App. A 2). These states are 
already eigenstates of and and a denotes a set 
of additional quantum numbers resulting from the cou- 
phng of the single spins s{i) to the total spin S. Due to 
Eqs. (Is]) the matrix elements of the Heisenberg Hamil- 
tonian (a' 5' M'| ^Heisenberg -5 M) between states with 
different S and M vanish, leading to a block-factorized 
Hamilton matrix in which each block can be diagonalized 
separately. In this case the field dependence of the en- 
ergies induced by the Zeeman term in Eq. ^ can easily 
be added without further complicated calculations. This 
is because the z-direction can be chosen to point along 
the external field, so that the Zeeman term commutes 
with ^Heisenberg, §^ and S^. Then M still serves as a 
good quantum number, and the effect of the external 
field B = B ■ on eigenstates of ^Heisenberg results in 
a simple field dependence of the energy eigenvalues Ei 
according to 



Ei{B)=Ei^g^BBMi . 



(9) 



This way thermodynamic properties depending on the 
temperature T and the external magnetic field B can 
easily be calculated from the energy spectrum of the in- 
vestigated magnetic molecule once the energies Ei are 
known. 



B. Irreducible tensor operator method 

The determination of the matrix elements of the 
Heisenberg Hamiltonian can elegantly be achieved with 
the help of irreducible tensor operatorsl ^^^^^ * To this end, 
it is necessary to reformulate the spin vector-operators in 



terms of irreducible tensor operators and to subsequently 
use tensorial algebra. In this regard the underlying the- 
ory is clearly based on group as well as representation 
theory. At this point it would probably not make sense 
to introduce all group-theoretical tools which lead to a 
complete understanding of the technical implementations 
used in this work. Several textbooks provide deep knowl- 
edge about these topics and the authors would like to 
refer to those-f^^H^ Nevertheless, at least the origin of 
appearing concepts and formulations shall be explained. 
Some understanding of abstract group and representa- 
tion theory is assumed. 



1. Irreducible tensor operators 

An irreducible tensor operator T^^^ of rank k is de- 
fined by the transformation properties of its components 
q under a general coordinate rotation R according to 

D{R)Tf^D-\R) = Y,TfD'^^,{R) . (10) 



Here D^^) denotes the operator associated with the co- 
ordinate rotation R. The subscripts q as well as q' take 
the values — /c, — A: + 1, . . . , /c. By D^^^J (R) the matrix ele- 
ments of the so-called Wig ner rotation matrices 'D^^^R) 
are denoted (cf. App. Al). 

The case /c = in Eq. (10) directly leads to what is 



called a scalar operator. As can easily be seen, a scalar 
operator is invariant under coordinate rotation, i.e. with 



5(i?)T(°)D-i(i?) = r(o)z)(°)(i?) 



(11) 



In analogy to the states which span the irreducible repre- 
sentation D^^^ of Rs and which are said to behave under 
coordinate rotations like the components of a vector the 
irreducible tensor operator T*^^^ is called a vector oper- 
ator. For example, the components of a first-rank irre- 
ducible tensor operator s*^^^ derived from the Cartesian 
components of the spin vector operator are given by 



.(1) 



(12) 



Stressing the analogy between the behavior of states 
and irreducible tensor operators under coordinate rota- 

(k) \ Iv 

tions, the role of the components Tq ^ in Eq. (10) has to 
be specified. The components of the irreduciEIe tensor 
operator T^^^ of rank k serve as a basis and therefore 
span the {2k + l)-dimensional irreducible representation 
D'^^^ of the rotation group Rs. 
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For comparison, consider a group Q and its irreducible 
representations r(^). The direct product of the irre- 
ducible representations F^*^ and F^-^^ separately spanned 
by two sets of basis vect ors is given by F^*^ F^-^^. It 
is reducible (cf. Eq. ( |A2| )) if linear combinations of the 
product functions can be found which transform as basis 
functions for an irreducible representation. This concept 
can - in a one-to-one correspondence - be extended to 
tensor operators since they behave like the above men- 
tioned functions. As a result, the direct product of two ir- 
reducible tensor operators spanning and 

can 

be decomposed into irreducible representations spanned 
by linear combinations of the products T^q^^^T^q^^\ The 
coefficients of these linear combinations are the Clebsch- 



Gordan coefficients of Eq. ( A4 ) . 

Formally, the direct proc 
operators is given by 



uct of two irreducible tensor 



{k) 



^qiq2q £qi £q2 ' V-L'^J 

<7i,<72 



where possible values of the resulting rank k can be de- 
termined in analogy to the vectorial coupling of spins 
and are given by /c = — /C2I, — fel + 1, • • • , + ^2- 
Equation ([l3| is a fundamental expression for the appli- 
cation of the irreducible tensor operator method within 
a numerical exact diagonalization routine. It leads to the 
desired formulation of the spin Hamiltonian in terms of 
irreducible tensor operators. 

As an example, the coupling of the first rank irre- 
ducible tensor operators U*^^^ and V*^^^ according to 
Eq. (13) shall be presented here. Considering a com- 
pound irreducible tensor operator with /c = 0, the cou- 
pling results in 



(0) 



V3 



(14) 



where for the Clebsch-Gordan coefficients the equation 
^liq% = VV^- {-^y~'^'^quq2 was uscd.^^ Exprcssiug 
the spherical components of V^^^ and U^^^ in terms of 



the Cartesian components in analogy to Eq. (12) yields 



(0) 



1 ^jjxyx _^ ^yyv _^ IfZyz^S 



(15) 



V3 



which is, apart from the prefactor, the scalar product of 
the Cartesian vector operators U and V. 

Finally, the problem of coupling irreducible tensor op- 
erators is identical to the coupling of angular momenta 
(cf. App. Al). Thus, from a mathematical point of 



view an advantage when using irreducible tensor opera- 
tors is that one can adapt the mathematical approaches 
for coupling angular momenta and that one can refer to 
them. 



2. Matrix elements of irreducible tensor operators 

In the case of an irreducible tensor operator T^^^ the 
matrix elements of this operator can be calculated ac- 
cording to the Wigner-Eckart theorem. It states for ma- 
trix elements with respect to spin states of the form 
\aSM) that 



{aSM\T[^^\a' S' M') = 
{-lf-^{aS\\T^^^\\a' S') 



S k (16) 
-M q M') ' 



The matrix element is apart from a phase factor 
decoupled into a Wigner-3J symbol and a quantity 
{a S\\T^^^ I \a' S') called the reduced matrix element of the 



irreducible tensor operator T^^^. The proof of Eq. (16) 



is given in standard textbooks about group theory and 
quantum mechanics. However, the physical meaning 
of this theorem and the consequences for the use within 
the irreducible tensor operator method shall be briefly 
discussed here. 

First of all, it should be mentioned that the Wigner- 
Eckart theorem relies on the transformation properties 
of the appearing wave functions and operators. Fur- 
thermore, since the reduced matrix element is com- 
pletely independent of any magnetic quantum number, 
the Wigner-Eckart theorem separates the physical part of 
the matrix element - the reduced matrix element - from 
the purely geometrical part reflected by the Wigner-3J 
symbol. The value of the reduced matrix element de- 
pends on the particul ar fo rm of the tensor operator and 
the states (cf. App. A 4) whereas the Wigner-3J sym- 



bol only depends on rotational symmetry properties. If a 
zero-rank {k = 0) irreducible tensor operator is assumed, 
the Wigner-3J symbol in Eq. (16) directly reflects that 



there is no transition between states \a S M) with differ- 
ent or M. The matrix of the Heisenberg Hamiltonian 
from Eq. (30), which can be written as a zero-rank irre- 



ducible tensor operator (see App. A3), therefore takes 
block-diagonal form without further calculations. 

According to Eq. ( 16 ) the calculation of the matrix el- 



ement of an irreducible tensor operator is directly related 
to the calculation of the reduced matrix element of this 
operator. The reduced matrix element of an irreducible 
tensor operator s^^^ with /c = 0, 1, acting on a basis func- 
tion of a single spin, can be derived from the evaluation 
of the Wigner-Eckart theorem. It yields the expressions 



(s||s(°)||s) = {s\\l}\s) = (2s + l)i 
(s||sW||s) = [s{s + l){2s + l)y^ 



(17) 
(18) 



where the zero-rank irreducible tensor operator of a sin- 
gle spin s*^^^ is given by the unity operator 1 and the 



components of s^^^ are given by Eq. (12). 

Using Eq. ([T3| in combination with tHe Wigner-Eckart 
theorem in Eq. (16) and a decomposition of states 
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\aSM) into product states according to Eq. (A4), one 
obtains the expression^^ 

{aiSia2S2S\\ jTC'i) ^tC'^)} ||a'i s'l 4 4 S") = 

(si s'l ki 
S2 S'2 k2 
S S' k 



(aiSi||T('=^)|Ks;)(a2S2||T('=^^IK4) • 



(19) 



This is the reduced matrix element of a compound irre- 
ducible tensor operator of rank k which consists of the 
direct product of two irreducible tensor operators of gen- 
eral ranks hi and /c2. 



Equation (19) is the basic formula for calculating re- 



duced matrix elements of irreducible tensor operators 
composed of single-spin tensor operators as they appear 
in the Heisenberg Hamiltonian. By a successive appli- 
cation every irreducible tensor operator of that kind can 
be decoupled into a series of phase factors, Wigner-9J 
symbols and the reduced matrix elements of single-spin 
tensor operators (cf. Eq. ([l2|). The successive applica- 
tion of Eq. (19) is often called decoupling procedure since 



the compound tensor operator that describes the system 
under consideration is decoupled so that its re duced ma- 
trix element can be calculated (see App. A 4). 



C. 



Point-group symmetries in Heisenberg spin 
systems 



Magnetic molecules, for instance those of Archimedian 
^ypglMmll are often - not only from a scientific point of 
view - perceived to be of a special beauty (cf. Fig. [5]). 
Certainly, this view is closely related to the high symme- 
try which can be found in the chemical structures and is 
referred to as point-group symmetry. Point-group sym- 
metries do not only contribute to the beauty of mag- 
netic molecules, but they are also very instrumental in 
characterizing the energy levels of the spectra and thus 
in extracting physical information from the underlying 
spin system. Often a numerical exact diagonalization re- 
mains impossible unless point-group symmetries are used 
in order to reduce the dimensionalities of the Hamilton 
matrices. 

It must be emphasized that there is a clear difference 
between a group of symmetry operations in real space 
that map the molecule on itself and the corresponding 
group of symmetry operations in a many-body spin sys- 
tem. Since from a physical point of view a magnetic 
molecule is described by a system of interacting spins 
with a certain coupling graph (cf. r.h.s of Fig. [T]), the 
term point-group refers rather to a group of symmetry 
operations on this coupling graph than to operations 
in real space. As a direct result, the group-theoretical 
characterization is also rather based on the symmetry of 
the coupling graph than on the molecular symmetry. Of 





FIG. 2: Sketch of the {Mo72Fe3o} molecule (l.h.s., picture 
taken from Ref . 60 ) . The underlying spin system exhibits the 
structure of an icosidodecahedron, i.e. possesses icosahedral 
{Ih) point-group symmetry (r.h.s.). 



course, since the number of appearing coupling constants 
as well as their strengths are estimated from the chemi- 
cal structure of the molecule, there is a close connection 
between the symmetries of the molecule and the corre- 
sponding coupling graph but not necessarily a one-to-one 
correspondence. 

In Heisenberg spin systems point-group symmetries 
can be included by mapping the symmetry operations 
on spin permutations. In order to emphasize this, 
the term spin-permutational symmetry instead of point- 
group symmetry is often used. However, in the work at 
hand the term point-group symmetry is used although 
the symmetry is always incorporated by mapping the 
point-group operations on spin permutations. In context 
to this, it has to be mentioned that in systems which in- 
clude anisotropics the incorporation of point-group sym- 
metries is much more complicated since rotations in real 
space have to be performed. 



CP « 
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FIG. 3: Coupling graph of a square of identical spins s with 
Da symmetry axes (l.h.s) and a rectangle with D2 symmetry 
axes (r.h.s.). 

Figure [3] exemplarily shows the coupling graphs of a 
square and a rectangle consisting of identical spins s with 
the symmetry axes of D/^ and D2 being indicated. The 
rectangle can be seen as resulting from the square - a 
system of four interacting spins with only one coupling 
constant J - by introducing a second coupling constant 
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of a diflFerent strength J' between spin pairs < 1 , 4 > and 
< 2,3 >. In this case the couphng strength is indicated 
by the length of the couphng paths between the spins. 
The introduction of the second couphng constant then 
results in a reduction of the point-group symmetry from 
to D2. 

The point-group operations on the spin system can be 
identified with permutations of the spins that leave the 
Hamiltonian invariant .^3 Such a permutation is repre- 
sented by an operator G{Ri) of the point-group Q for 
which the commutation relation 



H 



Heisenberg 



G{Ri) 







(20) 



holds, where i = 1, . . . , numbers the symmetry opera- 
tions up to the order h of Q. 

The theory of group representations now provides the 
theoretical background for the use of point-group sym- 
metries. The irreducible representations F^"^^ of a point- 
group G can be used to classify the energy eigenstates 
of ^Heisenberg and to block-factorizc the Hamilton ma- 
trices. The dimensionalities of the resulting subspaces, 
i.e. blocks, can be calculated with only little informa- 
tion. The irreducible matrix representations T'^^^Ri) of 
the group elements, i.e. the permutation operators, are 
somewhat arbitrary concerning the choice of the underly- 
ing basis. Thus, these elements are represented by their 
character, i.e. the trace of the particular matrix represen- 
tation. The character x^^\Ri) is invariant under unitary 
transformations and is in general given by 



{R) =TtT^''\R) = J2t^''\R) 



XX 



(21) 



A=l 



where In denotes the dimension of the n-th irreducible 
representation T^'^\ 

A given character table of Q now enables one to cal- 
culate the dimensions of the resulting blocks within the 
Hamilton matrix, i.e. the dimensions of the subspaces 
n{S, M, r(^)). Character tables for various point-groups 
can be found in almost every textbook about the the- 
ory of group representations. The authors would like to 
refer to these concerning the construction of character 
tables.l^^l^ 

In order to calculate the dimensions of the sub- 
spaces M, r*^^^), the reducible matrix representa- 
tions G{Ri) of the operators G{Ri) have to be considered. 
With /, m = 1, . . . , dimH(6', M) the matrix elements of 
these matrices are given by 

G{R^)im = {aiSM\G{R,)\amSM) , 

where the subscripts attached to a indicate that specific 
basis states are considered. The decomposition of the 
character xi^i) = G{Ri)ii with respect to the irre- 
ducible representations n of Q then yields 



(22) 



Using the great orthogonality theorem 
equation results in the expression 



dimH(5, M,r^^)) 



the above 



(23) 



where gives the multiplicity of the irreducible repre- 
sentation r^"^^ that is contained in the reducible repre- 
sentation G(C/c). C refers to the classes which the group 
elements can be divided in. Each class contains equiva- 
lent operations, for example n-fold rotations, which are 
linked by the same group operation and thus hold the 
same character. denotes the number of elements of 

From Eq. ([23| the dimensions of the subspaces 
can be calculated, but no information 
about the basis states that span these Hilbert spaces is 
given. The required symmetrized basis states can be de- 
termined by the application of the projection operatoif^ 



h ^ 



VQiR.) 



(24) 



This operator projects out that part of a function 
which belongs to the A^-th row of the irreducible represen- 
tation r*^^^. A subsequent application of the operator in 
Eq. on basis states, for example of the form la 5* M), 
is called the basis-function generating machine)^ 

Although Eq. ( 24 ) provides the information required to 



construct symmetrized basis states that serve as a basis 
for the irreducible representations of ^, it is important to 
notice that the matrices of F^^^ have to be known com- 
pletely. Of course, this does not cause a problem as long 
as Q entirely consists of one-dimensional irreducible rep- 
resentations. However, if more-dimensional irreducible 
representations appear, it is often more convenient to use 
the projection operator 



E^«"^ = iE[x^"H^^)rG(i?i) 



(25) 



that only requires information which can be extracted 
from the character table of Q. The operator of Eq. (25) 



projects out that part of a function \a S M) which be- 
longs to the irreducible representation T^'^\ irrespective 
of the row. As a consequence one has to orthonormal- 
ize resulting functions - for example by a Gram- Schmidt 
orthonormalization - in order to obtain the correct sym- 
metrized basis functions (cf. App. A 5). 



D. 



Point-group symmetry operations acting on 
vector-coupling states 



In the previous section II C some general remarks 



about the use of point-group symmetries have already 
been made. In order to transform the Hamilton matrix 
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to a block- diagonal form with respect to irreducible rep- 
resentations r*^"^^ of a point group one has to construct 
symmetrized basis functions \a S MT^'^'^). Equation (25) 
provides the projection operator that projects out that 
part of a state which belongs to the n-th irreducible rep- 
resentation of Q. However, the main challenge when con- 
structing symmetrized basis states is to find an expres- 
sion for the action of the operators G{Ri) on states of 
the form \aSM). 

The above mentioned operators G{Ri) which corre- 
spond to operations on the coupling graph (cf. Fig. [3| 
can be defined by their action on the product basis com- 
posed of the single-spin eigenstates of s^{i)' The states 
of the product basis of N identical spins can be denoted 
by 

l^imi) (g) |52m2) (8) • • • (8) \sNmN) \mi m2 . . . tun) 

and fulfill the eigenvalue equation 

s^{i) \mi . . . rrii . . . ttin) = rrii \mi . . . rrii ... itin) 

according to their definition. 

Now, in a spin system consisting of N spins an op- 
erator G{R) is considered that is given in the form 
G(7r(l)7r(2) ... 7r{N)) and describes a point-group op- 
eration. The notation G{7t) indicates for alH = 1, . . . , TV 
that the point-group operation interchanges the spin at 
site i with the spin at site 7r{i). The action of the opera- 
tor G{R) on the product basis is given by 

G{R) \mi 1712 ... rriN) \m^{i) m^{2) • • • ^7r(Ar)). (26) 



schemes a and h which is induced by an operator G rep- 
resenting a point-group operation: 



In Eq. (26) the following happens: by the action of G{R) 
the single^spin system at site i takes the z-component 
which the system at site 7r(i) has taken previously. The 
result is a different state of the product basis with per- 
muted m- values. 

The operators G{R) are only defined by their action 
on product states whereas the details of constructing 
symmetrized basis states that are linear combinations 
of vector-coupling states of the form \aSM) are still 
unknown. How to find an expression for the action of 
G{R) on these vector-coupling states is - to some extent 
- shown in Refs. [34] and [37l 

However, a slightly different perspective concerned 
with the application of general symmetry operations to 
a vector-coupling basis shall be presented here (for fur- 
ther technical details see App. A 5). In order to clarify 



the action of an operator representing a point-group op- 
eration on a vector-coupling state, the states are labeled 
according to the coupling scheme they are belonging to. 
Additionally, the set of quantum numbers referring to the 
coupling scheme is abbreviated by Greek letters. Primes 
indicate different basis states within the same coupling 
scheme. With these conventions one obtains the follow- 
ing general result for a transition between two coupling 



G \aSM), = 5^,p \pSM)k . 



(27) 



The Kronecker symbol 5cx^j3 on the right hand side indi- 
cates that the values of the spin quantum numbers of the 
different sets a and {3 are the same. Re-expressing the 
right hand side of Eq. (27) within states belonging to the 



coupling scheme a, i.e. inserting a suitable form of the 
identity operator j^, directly leads to the very important 
final result 

G \aSM)a = 

^^,,;3 a{aSM\pSM), \a'SM)a ^^^^ 

ol' 

which contains so-called general recoupling coefficients 
a{a S M\P S M)i). A general recoupling coefficient can be 
seen as a scalar product between vector-coupling states 
belonging to different coupling schemes a and b. 

It is by no means trivial to find an expression for a 
recoupling coefficient relating different coupling schemes 
if more than three or four spins are present. By defi- 
nition the expressions for the elements of the transition 
matrix relating two different coupling schemes result in 
Wigner-nJ symbols (cf. App. Al). While for three in- 



teracting spins Wigner-6J symbols occur, the size n of 
these Wigner coefficients increases with every additional 
spin taken into account. For four interacting spins the re- 
coupling coefficient is expressed by a Wigner-9J symbol 
and for five interacting spins the recoupling is described 
by Wigner- 12 J symbols. 

From a computational point of view it turns out to be 
unfavorable to use algebraic expressions for higher sym- 
bols than Wigner-9J symbols although there exist ex- 
pressions for 12 J symbols and 15 J symbols. In order to 
find an effective way to describe general recoupling coeffi- 
cients, one should use expressions in which only Wigner- 
6J symbols appear. These 6J symbols can be calculated 
using the formulsP^ 



A(ji, j2, j3)A(ji, J2, Js) 



A(Ji,j2, J3)A(Ji, J2,j3)^ 



fit) 



(29) 



where the triangle coefficient A(a,6, c) reads 

^ , /(a + 6-c)!(a-6 + c)!(-a + 6 + c)!) 

A(a,o, c) = ; -— 

^ ^ V (a + 6 + c+l)! 

The function f{t) in Eq. ([29]) is given by 



fit) = it- ji - 32 - h)\ it - 31 -J2- JsV- 

X{t-Ji- j2 - J3)! it-Jl-J2- jsV- 

X (jl + j2 + Jl + J2 - t)\ (j2 + j3 + J2 + J3 - ty. 

X ijs^ji^Js^Ji-ty. . 
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The sum in Eq. (29) is running over nonnegative integer 



values of t for which no factorial in f{t) has a negative ar- 
gument. Since even the evaluation of Wigner-6J symbols 
is a rather involved task as can be seen from Eq. (29), 



it is helpful to analyze the symmetry properties of the 
appearing symbols (cf. Fig. 14) in order to reduce the 
computational effort when constructing symmetrized ba- 
sis states. In this regard, when expressing the action 
of a point-group operation on a basis state according to 
Eq. ([28]), only those recoupling coefficients have to be 
calculated which are non-zero. 

So far, nothing has been said about the generation of 
a formula that describes a general recoupling coefficient. 
Finding a formula which only contains Wigner-6J sym- 
bols and some phase factors can most efficiently be done 
by using graph theory (see App. |B|. However, the choice 
of the coupling scheme which underlies the construction 
process of the basis states is essential for an effective 
computational realization of the numerical exact diag- 
onalization. In this regard the invariance of the coupling 
scheme under the applied point-group operations is a de- 
sired property (see Apps. |A 5| and A 6). 



III. APPLICATIONS 

In this section we like to present three applications for 
realistic spin systems of unprecedented size that can be 
treated using irreducible tensor operator techniques and 
point-group symmetries: a cuboctahedron of N = 12 
spins of spin quantum number 5 = 3/2 (Hilbert space 
dimension 16,777,216), an icosahedron of = 12 spins of 
spin quantum number 5 = 3/2 (same dimension), and a 
spin ring of A = 10 and s = 5/2 (Hilbert space dimension 
60,466,176) known as the ferric wheel Feio!^ 





FIG. 4: Structure of the cuboctahedron (left) and the icosa- 
hedron (right). 

It is very important to note that for the case of an- 
tiferromagnetic coupling the cuboctahedron as well as 
the icosahedron, see Fig. [4j belong to the class of geo- 
metrically frustrated spin systems^^ thus 
hardly accessible by means of Quantum Monte Carlo 
(QMC) calculations due to a so-called negative-sign 
problem .'^^HzH Complete diagonalization is therefore the 
only way to study the interesting behavior of such sys- 
tems both as function of field and temperature.^^ 44 | 62 | 77 | 



Density Matrix Renormalization Group (DMRG) tech- 
niques provide a very powerful approximatio n mainly for 
one-dimensional spin systems such as chains.CnilSlMThe 
method delivers the relative ground states for orthogo- 
nal subspaces. Extensions to include the approximate 
evaluation of excitations have been developed recently,!^ 
but are still primarily useful for one-dimensional systems. 
For the example of Feio, this method could provide low- 
est levels for each total magnetic quantum number M 
and thus one could evaluate the crossing fields of the 
lowest levels. More or less the same holds true for the 
Lanczos diagonalization technique, with which again 
extremal energy eigenvalues and eigenstates can be eval- 
uated. Quantum Monte Carlo would be able to evaluate 
observables for Feio, since this ten-site system is non- 
frustrated. Nevertheless, QMC cannot determine higher 
lying energy levels which would be crucial for e.g. inelas- 
tic neutron scattering studies. 

Since the three investigated systems exhibit a highly 
symmetric structure they can be modeled with just one 
nearest neighbor interaction, therefore the Heisenberg 
Hamiltonian simplifies to 



^Heisenberg = X] S(^) * * 



(30) 



<^,J> 



The summation is running over pairs < z, j > of nearest 
neighbors of single-spin vector operators s at sites i and j 
counting each pair only once. For the following examples 
J < is assumed. 



A. The cuboctahedron 

The magnetism of antiferromagnetically coupled and 
geometrically frustrated spin systems is a fascinating 
subject due to the richness of phenomena that are 
observed .'^^'^ One of the most famous spin systems is the 
two-dimensional kagome lattice .^^''^ It is very interest- 
ing and from the point of theoretical modeling appealing 
that similar but zero- dime nsional spi n systems - in the 
form of magnetic molecule j^^ * ^-*^*^^^ ^ - exist that poten- 
tially could show many of the special features of geomet- 
rically frustrated antiferromagnets. The cuboctahedron 
which consists of squares surrounded by triangles serves 
as one zero-dimensional "little brother" of the kagome an- 
tiferromagnet; the icosidodecahedron, which consists of 
pentagons surrounded by triangles, is another one. Such 
finite size antiferromagnets offer the possibility to dis- 
cover and understand properties that are shared by the 
infinitely extended lattices. An e xamp le is the discovery 
of localized independent magnons,'^^'^ which explain the 
unusual magnetization jump at the saturation field. Also 
the plateau at 1/3 of the saturation magnetization that 
appears in systems built of corner sharing triangles could 
be more deeply investigated by looking at the cubocta- 
hedron and the icosidodecahedron !^^*^^* 

The energy levels of the cuboctahedron {N = 12, 
s = 3/2, Hilbert space dimension 16,777,216), shown 
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FIG. 5: Full energy spectrum of a cuboctahedron with s = 
3/2. The energy levels have been calculated using D2 point- 
group symmetry. 



in Fig. [5) could be numerically evaluated using the D2 
point-group symmetry. For low-lying sectors of S = 
0,1,2,3 we also determined the energy levels according 
to irreducible representations of the full octahedral group 
Oh, see Fig. |6] One feature that can be clearly seen in 
Fig. |6] is the existence of an additional low-lying singlet 
below the first triplet. 



Ft 
□ 



n 



□ Alg 
A2g 
A Eg 

V Tig 
T2g 

Alu 
( A2u 
A Eu 

V Tlu 
<^ T2u 



the regular cuboctahedron with s = 3/2. This curve 
shows the aforementioned plateau at 1/3 of the satura- 
tion magnetization and a jump to saturation of height 
AM = 2. 
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FIG. 8: Zero- field heat capacity (top) and zero-field suscepti- 
bility (bottom) for the regular cuboctahedron with s = 1/2, 
s = 1, and s — 3/2. 



"^^0123 

s 

FIG. 6: Low- lying part of the energy spectrum of a cubocta- 
hedron s = 3/2. The energy levels are labeled according to 
irreducible representations of the full octahedral group Oh- 
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FIG. 7: Magnetization as a function of applied field at T = 
for the regular cuboctahedron with s — 3/2. The extended 
plateau at 1/3 of the saturation magnetization is clearly vis- 
ible. 

Figure [7] shows the magnetization curve at T = for 



Figure [s] compares the heat capacity (top) and the 
zero-field susceptibility (bottom) for the regular cuboc- 
tahedron with 5 = 1/2, s = l, and s = 3/2. The heat 
capacity shows a pronounced double peak structure for 
5 = 1/2 and s = 1 which dissolves into a broad peak 
with increasing spin quantum number. The broad peak 
also moves to higher temperatures with increasing s. The 
reason for the first sharp peak is twofold. Since there are 
several gaps between the low-lying levels the density of 
states has a very discontinuous structure which results 
in the double peak structure. For 5 = 1/2 the low- lying 
singlets provide a very low-lying non-magnetic density 
of states which is responsible for the fact that the first 
sharp peak is at such low temperatures. For 5 = 1 the 
first sharp peak results from both excited singlet as well 
as lowest triplet levels. For 5 = 3/2 a remnant of the 
first sharp peak is still visible; it is given by the low- lying 
singlets, but since they are so few, also influenced by the 
lowest triplet levels. The behavior of the heat capac- 
ity is contrasted by the susceptibility, bottom of Fig. [Sj 
which reflects mostly the density of states of magnetic 
levels and is only weakly influenced by low- lying singlets. 
Therefore, the first sharp peak, or any other structure at 
very low temperatures, is absent. 
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B. The icosahedron 

A spin system where the spins are mounted at the ver- 
tices of an icosahedron and interact antiferromagneticahy 
along the edges seems to be rather appeahng since it ex- 
hibits unusual frustrat ion p roperties such as metamag- 
netic phase transitions j ^^ * ^^" ^ Unfortunately, it appears 
to be challenging to synthesize such structures although 
icosahedra are rather stable cluster configurations!^ 




FIG. 9: Low-lying energy levels of an antiferromagnetically 
coupled icosahedron = 12 and s = 3/2. 

Here we investigate the properties of an icosahedron 
(A^ = 12) with single spin quantum number s = 3/2 
(Hilbert space dimension 16,777,216). The complete set 
of energy eigenvalues has been determined using D2 sym- 
metry. Figure [9] displays the low- lying part of all levels. 
Looking at the data file it turns out that very many (re- 
ally unusually many) levels are highly degenerate, which 
is in accordance with earlier investigations of ico sahedra 
with smaller single spin quantum number .^^^^^ We find 
that for every sector of total spin S and total magnetic 
quantum number M the irreducible representations A2, 
5i, and B2 contain exactly the same energy eigenvalues, 
whereas Ai is different. The non-degenerate ground state 
belongs to Ai. In addition very often near degeneracies 
occur. 



Figure [TO] shows the related zero-field heat capacity 
(top) and zero-field susceptibility (bottom). While the 
susceptibility does not look so unusual compared to other 
spin structures, the heat capacity appears to be really 
weired. Half way up the initial rise at very low temper- 
atures there is a Schottky-like peak that stems from the 
slightly split ground state {S = 0) levels. The further rise 
is due to the fact that the lowest 2x3 degenerate {S = 1) 
level is energetically rather close. Although higher-lying 
levels are separated by gaps from the lowest levels they 
also contribute to the heat capacity due to their massive 
degeneracy. Altogether this results in a low-temperature 
heat capacity that is much larger than the heat capacity 
of the cuboctahedron, compare Fig. [Sj 

The icosahedron may also serve as an example for the 
technical complexity due to the evaluation of recoupling 
coefficients. When combining 1^ symmetry with SU{2) 
many different recoupling formulas have to be generated 




FIG. 10: Zero-field heat capacity (top) and zero-field suscep- 
tibility (bottom) for the regular icosahedron with s = 3/2. 




FIG. 11: Low-lying part of the energy spectrum of an icosa- 
hedron with s = 1. The energy levels are labeled according 
to irreducible representations of the full icosahedral group Ih . 



(119 for the 120 group operations minus identity to be 
precise). This renders a treatment of the 5 = 3/2 icosa- 
hedron in Ih impossible. Although the sizes of the Hamil- 
ton matrices for the irreducible representations would be 
small, it is the construction of basis states that turns 
the evaluation of the needed matrix elements into a very 
lengthy procedure. Even in a parallelized job on 256 
cores on a supercomputer this task could not be com- 
pleted. For the 5 = 1 icosahedron we could finish a de- 
composition into irreducible representations of Ih, but - 
to give an example - the construction and subsequent 
diagonalization of the largest subspace, which has only a 
dimension of 3315, took about three days on 128 brand 
new Nehalem processors. Figure [TT] shows the low- lying 
part of the energy spectrum of an icosahedron with 5 = 1. 
The energy levels are labeled according to irreducible rep- 



11 



resentations of the full icosahedral group Ih- The exact 
and near degeneracies that have been discussed above 
for the case of s = 3/2 and D2 can now be resolved. 
For instance, the lowest S = 1 level belongs to T2u and 
is thus threefold degenerate. It is split from the higher- 
lying T2g by only 3/1000| J|, which in any calculation or 
measurement would look like a six-fold degeneracy, see 
also Ref. i95| 



C. Rings 

Molecular ferric wheels are among the very first mag- 
netic molecules. It appears that they can be synthe- 
sized in man y even nu mbered sizes, e.g. TV = 6, 8, . . . , 18 
Fe(III) spins.CSHHinil Since the spin of the Fe(III) ions 
is s = 5/2 the Hilbert space dimension grows rapidly 
from one ring size to the next. For the ferric wheel Feio 
it reaches already 60,466,176 rendering a complete di- 
agonalization impossible, at least in the past.^^ Based 
on the observation that the lowest energy eigenvalues 
Emin{S) obey a q uadratic dependence with respect to 
total spinpHIi3ll03j[i06] ^j^i^h is Lande's interval rule, ap- 
proximations could be derived for the level crossing 
fields^^^ as well as for low-lying excitations in for instance 
the truly giant ferric wheel FeigP^ As mentioned earlier, 
QMC is also capable of evaluating observables for even- 
membered unfrustrated spin rings. Nevertheless, none 
of the methods is able to determine higher-lying states. 
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FIG. 13: Susceptibility (top) and magnetization (bottom) of 
an antiferromagnetically coupled spin ring with = 10 and 
s = 5/2. The exchange parameter J = —9.6 cm~^ as well as 
the susceptibility data are taken from Ref. |70l 




FIG. 12: Low-lying energy levels of an antiferromagnetically 
coupled spin ring with = 10 and s = 5/2. 

In the following we present the first complete diago- 
nalization study of a spin ring similar to Feio, i-e. with 
= 10 and s = 5/2. This enables us to subsequently 
evaluate all thermodynamic functions, all excited levels, 
and if needed even the evolution for time-dependent ob- 
servables. For the diagonalization we used only the D2 
symmetry because it reduces the matrices already suf- 
ficiently and allows a faster computation of recoupling 
coefficients and thus matrix elements compared to the 
Cio symmetry. 



Figure 12 displays the low- lying energy levels for var- 
ious sectors of total spin S. The rotational band struc- 
ture of energy levels, which is at the heart of the afore- 



mentioned approximation, is clearly visible. Having ob- 
tained the eigenvalues an evaluation of the magnetization 
as function of temperature and field is easily possible. 
Figure [13] shows the susceptibility as function of temper- 
ature (top) and the magnetization as function of applied 
field (bottom) of an antiferromagnetically coupled spin 
ring with = 10 and s = 5/2. The exchange param- 
eter J = —9.6 cm~^ as well as the susceptibility data 
are taken from Ref. [70l ^ = 2.0 was assumed for the 
calculation as in the original paper. Since the theoreti- 
cal susceptibility curve is now known exactly, a detailed 
discussion of the experimental data becomes possible. As 
conjectured already in Ref.^^the low-temperature hump 
of the susceptibility is not a feature of the ring Hamil- 
tonian. We found numerically that it can also not be 
produced by a (reasonable) next-nearest neighbor inter- 
action. The authors of the experimental paper suggested 
that it might be produced by dimers of Fe(III), that 
should also contribute (slightly) to the large tempera- 
ture behavior, where the agreement with the theoretical 
curve is also not perfect. Our obtained energy eigenval- 
ues would allow now to set up a model that contains the 
ring and an unknown amount of dimer impurities. Open 
parameters would then be the amount of impurity and 
the g value which very likely deviates slightly from 2. 
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IV. OUTLOOK 

In this review we have undertaken the attempt to 
explain how numerical approaches to diagonalize the 
Heisenberg Hamiltonian work, that employ the spin- 
rotational symmetry SU{2) in combination with point- 
group symmetries. We hope that we could provide the 
reader with detailed knowledge on how to set up such a 
diagonalization scheme, especially through the extended 
technical appendix. 

What remains open for the future is to develop effi- 
cient schemes to evaluate recoupling coefficients, which 
at the present stage can be a very time consuming proce- 
dure that sometimes renders the calculation of matrix el- 
ements practically impossible although the complete ma- 
trix would be rather small in the end. A natural next step 
consists in finding optimal coupling schemes for a given 
molecule and point-group symmetry. 
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Appendix A: Realization of the irreducible tensor 
operator technique 

In this section the theoretical foundations presented 
in Sec. [ll] shall be specified. To this end, after having 
developed a basic idea of the coupling of angular mo- 
menta a spin square deals as a small example system. It 
is demonstrated how an appropriate basis can be con- 
structed and how the Heisenberg Hamiltonian looks like 
when expressed in terms of irreducible tensor operators. 
Further on, it is shown how the matrix elements can be 
evaluated by decoupling the irreducible tensor operator 
that describes the system. A central aspect is the use of 
point-group symmetries in combination with irreducible 
tensor operators that leads to the occurrence of general 
recoupling coefficients. 



1. Coupling of angular momenta and Wigner-nJ 
symbols 

As a first step the coupling of two general angular mo- 
menta j 1 and j 2 shall be discussed. Regarding this. 



the work at hand mainly refers to the use of defini- 
tions and name conventions that have been introduced 
by Wigner.l^ 

The vector coupling rule for the addition of angular 
momenta, known from elementary atomic physics, states 
that the resulting angular momenta J can be character- 
ized by a quantum number J. J assumes all values 

b'l - 32\Ml - j2| + + j2 - + j2 , (Al) 

where ji and j2 denote the quantum numbers of the an- 
gular momenta j i and j 2- Equation ( Al ) is also referred 

to as triangular condition for the coupling of two angular 
momenta. 

From a group-theoretical point of view the eigenstates 
\j m) of and j ^ span a (2j + l)-dimensional irreducible 

representation D^^^ of the group i?3, i.e. the group of all 
rotations within three-dimensional space. Keeping this in 
mind, the above vector addition rule Eq. (Al) is equiva- 
lent to 



J1+J2 

E 

j=\h-h\ 



(A2) 



According to this equation, the direct product between 
two irreducible representations of the dimensions (2ji + 
1) and (2^2 + 1), which is in general reducible, decays 
into (2 J+l)-dimensional irreducible representations with 
respect to J. 

Now, the operators D^^^\R) and D^^'^\R) are asso- 
ciated with an arbitrary coordinate rotation R of R^ 
and operate in the Hilbert spaces spanned by |jimi) 
and 1^2 ^2)- What does Eq. ( |A2| ) mean for the direct 
product T>^^^\R) T>^^'^\R) of the corresponding ma- 
trix representations? The product can be transformed 
into a block-diagonal form U(i?) by a unitary matrix A, 
i.e. 

'D^^'\R)^T>^^^\R) = AtU(i?)A , (A3) 
where U(i?) has the form 
V{R) = 

"D(l^'i-^2l)(i^) ... 

D^l^'i-^2l+i)(i?) ... 

... D(^'i+^2)(i?)_ 

The matrices D^^^(i^) appearing therein comprise matrix 
elements with respect to those functio ns th at span the 
irreducible representations D^^^ in Eq. (A2). 



The determination of the elements of the transforma- 
tion matrix A has been a central task in the theory of 
group representations. The elements of A are the so- 
called Clebsch-Gordan coefficients C^^^^^j^. They ap- 
pear in a more familiar form as scalar products between 
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a state of the form \ji j2 J M) and the product states 
I ji ^1 32 ^2) leading to the decomposition 

\j1j2JM) = {jimij2m2\jij2JM)\jimij2m2) 

mi ,7712 

= E C^^,'l^iM\jlm^j2m2) . (A4) 

rUi ,7712 

The Clebsch-Gordan coefficients therefore relate two dif- 
ferent orthonormal sets of basis vectors. It should be 
emphasized that these sets are obviously not orthogonal 
to each other because they span the same space. The 
Clebsch-Gordan coefficients are non-zero only if the vec- 
tor addition rule from Eq. ( |A1[ ) and simultaneously the 
equation mi + m2 = M hold. A very important symme- 
try of the Clebsch-Gordan coefficients is 



^7711 7772 M 



V ^) ^1712 mi M 



(A5) 



±1 according to Eq. (Al). 



with {-ly 

In order to reveal the symmetry properties of the 
Clebsch-Gordan coefficients, they are reformulated in a 
straightforward manner. A proper reformulation leads to 
the Wigner coefficients or Wigner-3J symbols 



ji h J 
mi m2 M 



which are related to the Clebsch-Gordan coefficients by 



ji h J 
mi m2 M 



(A6) 



-32 



-M 



The relation between Clebsch-Gordan coefficients and 



the Wigner-3J symbols from Eq. (A6) directly leads 



to non-zero values of the Wigner-3J symbols only if 
mi + m2 + M = holds and the vector addition rule 



from Eq. (Al) is fulfilled. 



Expressed in terms of Wigner-3J symbols the symme- 
try property of the Clebsch-Gordan coefficients given in 
Eq. ( A5 ) takes the form 



ji h J 
mi m2 M 



(^_iyi+32+J 



h ji J 

1712 ^1 M 



(A7) 



Thus, the Wigner-3J symbols are invariant under an even 
number of permutations of the columns whereas under a 
single permutation they obey Eq. (A7). A further evalu- 



ation of the symmetry properties of the Clebsch-Gordan 
coefficients yields an additional symmetry of the Wigner- 
3J symbols. 



I.e. 

ji 32 J 
mi m2 M 



(-1) 



3 1+32+ J 



Ji 

-mi 



J2 

-m2 



-M 



(A8) 



So far the coupling of only two angular momenta has 
been considered. A procedure similar to the one which 



has led to the Wigner-3J symbols now leads to the oc- 
curence of Wigner- 6 J symbols. In the case of a coupling 
of three angular momenta, J = ji + j2 + j35 a basis 

can be constructed in which the representations of the 
operators and as well as ji, and are diag- 
onal. Obviously, there exists a certain freedom of choice 
in the construction of this basis. The resulting J can be 
constructed in three different ways, namely 



J = ( j 1 + h) + js , 
J = Ji + (j2 + js) , 

J = (jl + js) + j2 . 



(A9) 
(AlO) 
(All) 



This leads to three different basis sets, each with the 
square of one of the operators j' = ji + j2, j^' = j2 + js, 

and = j 1 + j s given in a diagonal form. The ma- 
trix elements of the unitary transformation matrix which 
connects two of these sets of basis states can be found as 
scalar products between states belonging to two differ- 
ent coupling schemes. Expressing a state belonging to 
a coupling scheme resulting from a coupling according 
to Eq. (A9) in ter ms of states belonging to the scheme 
resulting from Eq. ( |A1Q ) yields 



\jl j2 jl2 js J M) = 

X^O'l j2 js j2S J\jl j2 jl2 js J) \jl j2 js j2S J M) , (^^2) 



323 



with the quantum numbers ji2 and j2s referring to j 
and y'^ ^ respectively. 



/2 



Scalar products of the kind found in Eq. (A12) - often 
called recoupling coefficients - are independent of any 
magnetic quantum number. They can be evaluated by 
decomposing the vector-coupling states into a sum of 
prod uct s tates with the help of an extended version of 
Eq. (A4). For example, the decomposition of the ket 



on the left hand side of the aforementioned recoupling 
coefficients yields 

b'l j2 js j2S JM) = 

^ O'l mi j2 m2 js ms \ji j2 js j2s J M) (^^;l3) 

mi ,7772 ,m3 

X \ji mij2m2jsms) . 



The scalar products in Eq. ( |A13 ), i.e. the matrix el- 
ements of the transformation matrix that connects the 
vector-coupling state with the product states, are called 
in analogy to the former name convention generalized 
Clebsch-Gordan coefficients Ct^^L^^l^^^, for the cou- 

''^ mi 77l2 m^ 77123 -'^^ 

pling of three angular momenta. They can be simplified 
to a product of Clebsch-Gordan coefficients according to 



/-ijl 32 33 323 J _ /-ij2 33 32 

^7771 7772 W3 77723 M / j ^7772 7773 



'-"777i7n23M • (A14) 
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As one can see, generating generalized Clebsch-Gordan 
coefficients for the coupling of more than three angu- 
lar momenta and - in addition - an expression of it in 
terms of Clebsch-Gordan coefficients is then a straight- 
forward task. Here it should be mentioned that accord- 



ing to Eq. (A14) generalized Clebsch-Gordan coefficients 



can also be expressed as a sum over products of Wigner- 
3 J symbols using the relation from Eq. (A6). Coeffi- 



cients of this kind are then called generalized Wigner 
coefficients—^ 



The recoupling coefficients that appear in Eq. (A12) 



can now be reformulated in terms of Wigner-6J symbols 
in order to reveal their symmetry properties. For the 
transition between the basis sets belonging to ji2 and 
j23 the corresponding Wigner-6J symbol is related to the 
recoupling coefficient hy^ 



\h J j23 J 



(-1) 



jl-\-j2-\-j3-\-J 



(2ji2 + l)- 



(A15) 



X (2j23 + 1) 2 {jl 32 js 323 J\3i 32 312 js J) • 



Equations (A13) and (A14) directly show that a Wigner- 
6J symbol can be expressed as a sum over products of 
Clebsch-Gordan coefficients or, by using Eq. (A6), over 
products of Wigner-3J symbols. 

For the sake of completeness, also the Wigner-9J 
symbol^^ shall be given which result as elements of the 
transition matrices when recoupling four angular mo- 
menta. For example the transition between two different 
sets of basis states yields a Wigner-9J symbol like 



Jl 32 312 
3s j4 jsA 

[(2ji2 + l)(2j34 + l)(2ji3 + l)(2j24 + 1)]- 
X O'l j2 312 js 34 3SA J\jl js jlS j2 ^4 ^24 J) 



(A16) 



A very comprehensive overview of algebraic expressions 
for Wigner-nJ symbols as well of their symmetry prop- 
erties is given in Ref. "58 , Regarding the use of Wigner 
symbols in connection with irreducible tensor operators, 
here only a graphical visualization of the triangular con- 
ditions for a Wigner-6J symbol is shown. The Wigner-6J 
symbol is nonzero only if for certain triads of quantum 
numbers the triangular condition (Eq. (Al)) holds. For 
which combination of quantum numbers of angular mo- 
menta the triangular condition has to be valid in order to 
yield a nonzero Wigner-6J symbol, is visualized graphi- 



cally in Fig. 14 

Since a Wigner-9J symbol can be written as a sum 
over products of Wigner-6J symbolJ^, these triangular 
conditions play an important role in reducing the com- 
putational effort when matrix elements of the Heisenberg 
Hamiltonian are calculated using irreducible tensor op- 
erators. 




FIG. 14: Graphical visualization of the four triangular condi- 
tions occurring in a Wigner-6J symbol which have to be valid 
in order to yield a nonzero result. 



2. The construction of basis states 

The reduction of the dimensionalities of the Hilbert 
spaces in which the Hamilton matrices are set up in or- 
der to solve the eigenvalue problem is always a desirable, 
but also numerically involved task. Especially if the ba- 
sis a priori reflects symmetry properties of the system, 
an appropriate choice can be of great help. As men- 
tioned in Sec. a basis that consists of vector-coupling 
states and incorporates full spin-rotational {SU{2)) sym- 
metry would be the ffist choice. In isotropic spin systems 
as described by a Heisenberg Hamiltonian the Hamilton 
matrix is then block-diagonalized with respect to S and 
M without further calculations since there are no tran- 
sition elements between states of a different total mag- 
netic quantum number M and different total-spin quan- 
tum number S. 

The vector-coupling states used in the present work 
appear in the form \aS M). a denotes a set of interme- 
diate quantum numbers resulting from the chosen cou- 
pling scheme according to which the spins are coupled. 
As mentioned in the previous section the choice of the 
coupling scheme is somewhat arbitrary since it only re- 
flects the bracketing in the expression for the total spin 
operator of the system S = s^. 

The simplest choice of a coupling scheme would prob- 
ably be a successive coupling of the single-spin vector 
operators s^. In the case of a spin square the set of in- 
termediate quantum numbers, if coupled according to a 
successive coupling scheme, looks like 

{Sl2 Si,Si23 S2} 

leading to a vector-coupling state of the form 
\s1S2S1SsS2S4SM). Here, the notation of the intermedi- 
ate spin quantum numbers, i.e. 512 and 5123, is changed 
in comparison to App. |A 1[ Intermediate spins are now 
numbered with respect to their order of appearance in 
the coupling scheme and additionally overlined. This no- 
tation has a clear advantage if larger spin systems are 
investigated and is used in the following. In order to 
clarify which spins are coupled, the single-spin quantum 
numbers Sj can also be found in the ket. It would not 
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be necessary to include them since they appear as fixed 
numbers, but it turns out to be more convenient. 
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FIG. 15: Successive coupling of four spins s = 1. The (red) 
subscripts denote the multiplicity, i.e. the number of paths 
leading to the spin quantum number. The gray numbers refer 
to the spin quantum numbers of the coupled single spins. On 
the left the coupling scheme is indicated. 



Now, if the coupling scheme is chosen and the frame- 
work of the resulting basis states is fixed, one has to 
construct the basis states by finding those values o f th e 
intermediate spins that are valid according to Eq. (Al). 



This procedure can be visualized by constructing a cou- 
pling pyramid as it is shown in Fig. [15] In Fig. 15 
four spins with 5 = 1 are successively coupled in order 
to yield the values of the total-spin quantum number 
S. The (red) subscripts next to the quantum numbers 
of the intermediate spins denote the number of different 
paths leading to the same value for an intermediate spin 
quantum number, i.e. the multiplicity. The small (gray) 
numbers interrupting the lines connecting different in- 
termediate spin quantum numbers indicate the quantum 
numbers of the single spins, i.e. Si = 1. For the sake 



of clarity, on the left of Fig. [15] the underlying coupling 
scheme is given once more. 
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FIG. 16: Pairwise coupling of four spins s = 1. The (red) 
subscripts denote the multiplicity, i.e. the number of paths 
leading to the spin quantum number. The gray numbers refer 
to the spin quantum numbers of the coupled single or inter- 
mediate spins. On the left the coupling scheme is indicated. 

Of course, a successive coupling scheme is not the only 



possible way to couple the single spins to a total spin S. 
For example, the construction of quasi-classical states as 
described in Refs. [38l and [39l requires a different coupling 
scheme. There, a coupling scheme is chosen in which 
spins belonging to a certain sublattice are coupled in or- 
der to get the total sublattice spin that is afterwards 
coupled to the total spin of the system. In Fig. [16] a cou- 
pling pyramid for a different coupling scheme of a spin 
square is shown. The spins Si, S2 and S3, S4 are coupled 

to yield intermediate spins and S^, respectively. The 
intermediate spins are then coupled to the total spin S. 

The resulting multiplicity of states with the same 
quantum number S is obviously independent of the cho- 
sen coupling scheme. At this point it is important to real- 
ize once again that these states, although resulting from 
different coupling schemes, form basis sets of the same 
Hilbert space H. In the case of a square the coupling of 
four single spins 5 = 1 discussed above results in basis 
states that span the subspaces 1-L{S) with S = 0,J^. . . , 4. 



Writing this result in a different way using Eq. (22), one 
obtains for the direct product of the irreducible represen- 
tations D^^"^ of the single spins the following expression: 

(g) (g) D^^^ (g) D^^^ = 

3 • D^^^ + 6 • D^'^ + 6 • D^^^ + 3 • D^^^ + 1 • D^^^ . 

The knowledge of the dimensions of the resulting irre- 
ducible representations, i.e. subspaces 1-L{S)^ is a central 
task whenever performing numerical exact diagonaliza- 
tion. Deduced from a successive coupling of the spins, 
the dimensions d\ml-L{S) can be calculated by a simple 
recurrence formula. 

The number of paths leading to a certain combination 
(5, n) is denoted hy ds- n refers to the number of partic- 
ipating spins in each step and runs from 1 to A^. n can 
also be seen as labeling the rows of the coupling pyramid 
for a successive coupling (see Fig. 15). In the case of 



a homonuclear system with A^ spins 5, the multiplicity 
ds{S^n ^ 1) is given by 



min(S'+s,n-s) 

ds(5,n+l)= ds{S',n), 

S' = \S-s\ 



(A17) 



where S lies in the interval that is bounded by 

/ . .^ . ^ ^ / if 25(n + 1) even . . 

(ri^l)s>S>^r if25(n + l)odd ^^^^^ 



according to a vector coupling rule (cf. Eq. ( |A1| )). The 
number of paths leading to a certain quantum number S 
for one spin is just 



-{I 



(A19) 



The dimension of the Hilbert space 7^(5') is then given 
by dimniS) = dsiS, N) ■ {2S + 1). 
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The dimensions within a heteronuclear spin system, 
i.e. a system with different values of single-spin quantum 
numbers, are obtained along the same route. However, 
the range of valid values for S has to be calculated for 
each st ep of recurrence separately, in contra st to the use 
of Eq. (A18). Furthermore, the sum in Eq. (A17) would 
run from 5" = — | to min(S' + s^, "^27=1 ^^i) where the 
spin quantum number of every single spin Si is individu- 
ally labeled by the index i. 



3. The Heisenberg Hamiltonian expressed as 
irreducible tensor operator 

In order to determine the spectra of magnetic 
molecules as it is done throughout this work, the Heisen- 
berg Hamiltonian of Eq. ( 30 ) has to be expressed as irre- 



ducible tensor operator. In this section a general expres- 
sion for the Heisenberg Hamiltonian is presented which 
can be used as a starting point for the calculation of the 
energy spectrum using the irreducible tensor operator ap- 
proach. 

spin dimer - The first step in finding an expression for 
the Heisenberg Hamiltonian in the form of an irreducible 
tensor operator is done by considering a spin dimer. The 
Hamiltonian of the dimer takes the simple form 



^dimer = - J s(l) • s(2) . 



(A20) 



Now, using Eq. (13) the above Eq. (A20) can easily be re- 



produced. Since the Heisenberg term is given by a scalar 
product, the corresponding compound irreducible tensor 
operator is of rank k = 0. Using Eqs. (12) - (15) one 
finds the expression 

{s«(l)®sW(2)}^°^ 



\/3 



S(l)-S(2) ■ 



Thus, the tensorial form of the Heisenberg Hamiltonian 
of a spin dimer is 

^dimer = V3J {s(ini)®sW(2)y°^ . (A21) 

spin triangle - Since the tensorial form of i^dimer in 



Eq. (A21) describes a simple bilinear spin-spin interac- 



tion, an expression for a general Heisenberg Hamiltonian 
can now be developed. As a first extension a spin triangle 
is considered. The Hamiltonian that has to be converted 



to an irreducible tensor operator is 

^A=- J (s(l)-s(2) + s(2).s(3) 
+g(3) •£(!)) . 



(A22) 



In a very general form the successive coupling of three 
single-spin irreducible tensor operators of ranks fci, k2 
and fcs leads to 



s^^'\l) s(^^)(2)}^^'^ s(^3)(3)| 



ik) (A23) 



Here k denotes the rank of the resulting irreducible tensor 
operator and ki the rank of the intermediate (coupled) 
one. The ranks of the many-particle tensor operators are 
given by the coupling rules for spin quantum numbers 
known from the spin vector coupling. For example, the 
rank ki is given by /ci = — /C2I, — /C2I + 1, • • • , +^2 
with k being determined accordingly. It must be empha- 

(k) 

sized that T: includes all spin-spin interactions of a 
trimeric spin system, and thus it has to be specified in 
order to give the desired tensorial formulation of ^a- 

With s^^^ = 1 and the tensorial expression found for 



a bilinear coupling in Eq. (A21) one arrives at 



+T^.°^(0, 1,1,1)' 



= T^\{k^}dk^}\k^ = k, = 1)(.A24) 



<^,J> 



The notation of T''^\{ki}^ {ki}\ki = kj = 1) indicates 



that only the ranks of single-spin tensor operators i and j 
are chosen to equal 1 whereas the other tensor operators 
are of zero rank. The rank of the intermediate tensor 
operator /ci is fixed by the contributions of ki and k2 to 
the zero-rank tensor operator T^.^~^\ 

Following Eqs. (|A21[) and (|A24[), the Heisenberg 



Hamiltonian of a general spin system is given by 

H Heisenberg 

^ E J^JT^''H{k^}Ak^}\k^ = k, = 1) , 



(A25) 



where the irreducible tensor operator T^^^ ({/c^}, {/c^}) di- 
rectly depends on the investigated system and the chosen 
coupling scheme for the coupling of the single-spin tensor 
operators s^^^\i). 

As an example, the resulting expression for the gen- 
eral, tetrameric tensor operator Th in the case of a 
spin square shall be presented. The underlying coupling 
scheme of the single-spin tensor operators is chosen to 



17 





ki 


k2 


k3 


/C4 


ki 


k2 


k 


< 1,2 > 


1 


1 

















< 2,3 > 





1 


1 





1 


1 





< 3,4 > 








1 


1 











< 4,1 > 


1 








1 


1 


1 






TABLE I: Values of the ranks for a spin square described 
by a Heisenberg Hamiltonian. The rows refer to the nearest- 
neighbor interaction < i,j > between single spins i and j. 



be pairwise because it is often advantageous to couple 
this way when add itionally point-group symmetries are 
used (see App. A 5). Then, the general irreducible tensor 



operator of a tetrameric system takes the form 
i('=i)(l) (g) s('=^)(2)y 



0{s(^^)(3)0s(^^)(4)} 



(fei) 



(A26) 



The values of the ranks appearing in Eq. (A26) for the 
spin square modeled by a Heisenberg Hamiltonian are 
shown in Tab. [Tj The ranks of the single- spin tensor op- 
erators are fixed by the particular spin-spin interaction. 
Then, the ranks of the intermediate tensor operators ki 
and k2 can be constructed from the known coupling rules 
with the given ranks /c^, i = 1, . . . , 4, so as to yield an ir- 
reducible tensor operator with k = 0. 



4. Matrix elements - decoupling 



The calculations of the matrix elements of the Hamil- 



tonian in Eq. (A25) are performed with the help of the 
Wigner-Eckart theorem fMIII Por this reason, the cou- 
pling scheme of the basis states \a S M) cannot be cho- 
sen independently from T^^^. Those quantum numbers 
which appear within the set a should also appear in the 
sets {ki} and {/c^}, i.e. the couplings for the generation 
of the basis states and the general irreducible tensor op- 
erator of the system should be chosen to be equal. Other- 
wise, transformations between states of different coupling 
schemes would be necessary. 

In order to show how matrix elements can be calcu- 
lated, the application of the decoupling procedure is dis- 
cussed for a square system. The Heisenberg Hamiltonian 
in the tensorial form can be derived from Eq. ( A25 ) and 
is given by 



Ha 



V3J J2 T^°H{ki},{ki}\ki 



kj 



1) 



(A27) 



The general irreducible tensor operator of the tetrameric 
system T^:^^ based on a pairwise coupling scheme was 
already presented in Eq. ( |A26| ). The values of the ranks 
appearing in T^f^~^\{ki} ^ {ki}) are tabulated in Tab. I 
The pairwise coupling scheme in the used construction oi 
T^f ^ corresponds to basis states of the form 

\aSM)=^ \s1S2S1SsS4S2SM) . 

By the application of the Wigner-Eckart theorem the 
calculation of the matrix elements of is now - apart 
from the prefactor and the summation over t he nearest- 
neighbor interactions according to Eq. (A27) - reduced 
to the evaluation of terms like 



{aS M\T[^\{h},{k^})\a' S' M') = 
i-M I M') 



In the case of a square four such terms appear in 
differing from each other by the values of {ki} and {ki}. 



Since the Wigner-3J symbol in Eq. (A28) is reduced to 



SOS' 
-M M' 



V2S~- 



MM' 



Eq. ( |A28[ ) itself is reduced to 
{aSM\T^\{h},{ki})\a' SM) = 

(_l\2iS-M) _ _ _, _, 

^ ^ - • (5152^15354^2^1 |T(°) 1 15152^15354^2^) • 

(A29) 



V2S~- 



Since all matrix elements between states of different S 
and M vanish, it becomes obvious that all calculations 
can be performed in subspaces 1-L{S^M). Furthermore, 
because H\j is independent of an external magnetic field 
and thus no M-dependence of the energies is given, the 
spectrum can be evaluated in subspaces 1-L{S^ M = S). 



Following Eq. (A29) the matrix elements can directly 
be obtained by determining the reduced matrix elements 
of Tf:^\ By a successive application of Eq. (jioj), i.e. a 

successive decoupling of the reduced matrix ele- 

ments are traced back to the reduced matrix elements of 



single-spin tensor operators that are given in Eqs. (17) 
and ([18| for /c = 0,1. 

(k) 

For the reduced matrix element of Th the decoupling 
yields 



<i,j> 
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{siS2SiSsS^S2S\\T[^^\\siS2S[ssS^S'^S') 




(fcl) 



l)]^[(25i + l)(25; + l)(2^i + l)]^ 



(A30) 



In the case of the appearing ranks of the tensor 
operators can be found in Tab. |l] and is forced to be 
S' = Shy the Wigner-3J symbol in Eq. (|A28D. 



The clear structure of the resulting expression for a 
reduced matrix element in Eq. ( A3Q| ) now allows one to 
write a highly flexible and structured computer program 
that takes over the calculation and diagonalization of the 
Hamilton matrix. Regarding the calculation of the ma- 
trix elements, a recurrence formula can be implemented 
which decouples the irreducible tensor operator of the 
system step-by-step 



5. Using point-group symmetries 

General considerations concerning the use of point- 
group symmetries in Heisenberg spin systems have al- 



ready been presented in Sees. |II C| and |IID[ Now, as a 
clarification of these considerations, a 7r-rotation around 
the central (72-axis of a spin square is considered (cf. 
Fig. 17). According to a successive coupling scheme. 



the vector-coupling basis states are given in the form 
\s1S2S1SsS2S4SM) . It must be emphasized here that 
the underlying coupling scheme, that determines the way 
how basis states are constructed, can be chosen inde- 
pendently from any symmetry considerations, although 
a suitable choice will reduce the calculations as is dis- 
cussed below. 
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FIG. 17: Sketch of a spin square with D4-symmetry opera- 
tions. 



Figure [T7| shows the coupling graph of the square with 
i^4-symmetry operations. The operations are labeled 
with respect to n-fold rotations around the given axes. 





Classes 


E 


C2 


2C4 


2^ 




op's 


5(1234) 


g(3 412) 


5(412 3) 
5(2341) 


5(214 3) 
5(4321) 


5(1432) 
5(3 214) 



TABLE II: Classes of the D4 point group with corresponding 
operators for the spin square. 



Operations belonging to the same class are marked with 
identical labels while identical operations belonging to 
different classes can be distinguished by the number of 
primes. 

In Tab. [IT] the operators corresponding to the D4- 
symmetry operations on the square are given. The oper- 
ators are classified with respect to the classes they belong 
to. The mentioned operator that acts in spin space and 
corresponds to a 7r-rotation around the central axis takes 
the form G(3412). 

According to the generalized Clebsch-Gor dan c oefii- 
cients for the coupling of four spins and Eq. (A14), un- 



folding the state \s1S2S1SsS2S4SM) into a linear combi- 
nation of product states results in the expression 



\s1S2S1SsS2S4S M) = 



Si S2 Si 

mi 1712 Ml 



C- 



S1S3S; 



-yS2 S4 S 



3 'J2 . 

ms M2 M2 m4 M 



=M 



(A31) 



X |mi m2 ms m^) 



The summation indices are entirely determined by the 
constraint rrii = M. The values of the intermediate 
magnetic quantum numbers Mi can be deduced from the 
magnetic quantum numbers of the involved single spins, 
i.e. Ml = mi + m2 and M2 = Mi + 7713. 

Following Eqs. (26) and ( A31[ ), performing the tt- 
rotation described by the operator G(3 412) results in 
the expression 



5(3 412) \s1S2S1SsS2S4SM) 

-isi S2 Si r^Si S3 S2 



~lS2 S4 S 



mi m2 Ml Ml ms M2 M2 m^ M (A32) 

T.^^^=M 

X |m3 7714 mi 7x12) . 

Due to the performed permutation on the product states 
the resulting state cannot easily be represented as a 
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vector-coupling state belonging to the former coupling 
scheme given by the successive addition of spin oper- 
ators: s(l) + s(2) = S(l), S(l) + s(3) = S(2), and 

S(2) + s(4) = S. 

At this point, it becomes obvious that the operator 
G is inducing a transition from the coupling scheme, 
according to which the basis states have initially been 
constructed, to another one. A proper re-labeling of 
the summation indices in the sum of Eq. (A32) with 
respect to a point-group operation ^ i.e. a (— tt)- 
rotation around the central C2-axis, reveals the resulting 
coupling scheme in which G(^A12)\siS2S 18^828 4,8 M) 

can be represented as a vector-coupling state. In this 
special case one finds that G(34 1 2) is inducing a transi- 
tion to a coupling scheme given by s(3) + s(4) = S(l'), 
S(r) + s(l) = S(20, and §(2^ + s(2) = S. As a short- 
hand notation one can write 



\siS2SiS^82SaSM) 



G(3412) 



\s^S4^Sv S182' S28M) 



with the limitation that this expression does not give the 
concrete values of the appearing quantum numbers of the 
states. Nevertheless, it illustrates the transition between 
vector-coupling states of different and thus independent 
coupling schemes. 



Following Eq. (28) the action of the 7r-rotation on a 



vector-coupling state of the chosen form results in 
G(3412) \siS2Sis^S2SaSM) = 

Yl ^S^,S,.^S2,S,^^l^2S[ssS2S4SM) 
X {siS2S[ssS2S4SM\ssS4SrSiS2'S2SM) . 



(A33) 



As mentioned in Sec. Ill Dl the main task when calcu- 
lating the action of a point-group operation on a vector- 
coupling state is the determination of the general recou- 
pling coefficients connecting states of the initial and the 
resulting coupling scheme. Generating a formula for gen- 
eral recoupling coefficients can only be performed in a 
rather advanced procedure (see Sec. IbI). 
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FIG. 18: Visualization of two possible couplings for the 
square: successive coupling scheme (l.h.s.) and pairwise cou- 
pling scheme (r.h.s.). 

Regarding the recoupling coefficients that appear when 
performing point-group operations, a very helpful simpli- 
fication shall be mentioned here. From Eq. (A32) it can 



be seen that the action of the operator G(34 1 2) on the 
product states prevents the re-expression of the resulting 
linear combination of product states as a simple vector- 
coupling state belonging to the initial coupling scheme. 
However, the choice of the coupling scheme according to 
which the initial basis was constructed is somewhat ar- 
bitrary. In order to minimize the computational effort, 
which is directly related to the numb er of states with 
non-zero recoupling coefficient in Eq. (A33), one has to 

A favorable 
This 
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choose a non-successive coupling scheme, 
coupling scheme of this kind is shown in Fig. 
scheme is referred to as pairwise coupling scheme and the 
basis states look like \s1S2S1SsS4S2SM). The 7r-rotation 
around the central C2-axis of the square now induces a 
transition that can be symbolized by 

_ _ G(3412) _ _ 

\s1S2S1SsS4S2SM) — > \ssS4Si'SiS2S2'SM) 



and leads according to Eq. ( |27| ) to a recoupling coefficient 
of the form 

{siS2SiSsS4S2SM\ssS4^SrSiS2S2^SM) . 

Now, the calculation of a formula for this recoupling co- 
efficient is trivial since the intermediate spin operators of 
the initial and the resulting coupling scheme are mutually 
the same, i.e. S(r) = S(2) and §(2^ = S(l). Unfolding 

\ssS4{Si' = S2)siS2{S2' = Si)SM) into a linear combi- 
nation of product states in analogy to Eq. (A31) leads 
to 



\s3S4S2S1S2S1SM) 



E 



^S3 S4 S2_ 

ma 7714 M2 



S2 Si G^'^ ^ 

mim2Mi M2M1M 



^. nii=M 

X |mi m2 ms 1714) . 

In order to convert the Glebsch-Gordan coefficients to 
a form that appears when unfolding states of the form 
\s1S2S1SsS4S2SM) one simply has to use the symme- 
try property of the Glebsch-Gordan coefficients from 
Eq. (|A5|. This leads to 



G 



S2 Si s 
M2 Ml M 



(-1) 



Si+S2- 



-S ^Si S2 s 
Ml M2 M 



and thus an expression for the recoupling coefficient is 
obtained which only contains one simple phase factor: 



{siS2SiSsS4S2SM\ssS4SvSiS2S2'SM) = {-lfi+S2-S 

The action of the operator performing a 7r-rotation on 
a state belonging to the pairwise coupling scheme men- 
tioned above therefore directly results in a state belong- 
ing to the same coupling scheme with an attached phase 
factor. Thus, it has been shown that with a cleverly 
chosen coupling scheme the computational effort, that is 
required for the calculation of symmetrized basis states 
according to Eq. (25), can be minimized. The graphi- 



cal visualization of possible coupling schemes in a square 
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shown in Fig. 18 makes it obvious that one wih find 
a simple recouphng formula depending only on phase 
factors whenever one can find a coupling scheme that 
is invariant under the performed symmetry operation. 
Since one has to sum over all symmetry operations of the 
group in order to construct symmetrized basis states of a 
given irreducible representation (see Eq. ([25|)), the under- 
lying coupling scheme should be chosen in such a way as 
to simplify all or at least most of the resulting recoupling 
coefficients. This means that the coupling scheme should 
be invariant under all or at least most of the symmetry 
operations. 

However, one will not always be able to find a coupling 
scheme that simplifies the calculation of the recoupling 
coefficients as shown above. Especially if the system un- 
der consideration is exhibiting three- fold symmetry axes, 
such a procedure turns out to be impossible by means of 
a pairwise coupling scheme. In this case a generalization 
of finding a recoupling formula independent of the choice 
of the coupling scheme becomes necessary. 



6. Computational effects of the choice of the 
coupling scheme 

The computational realization of the theoretical back- 
ground presented in this work has been a central task. 
The performed calculations would not have been pos- 
sible without developing a highly parallelized computer 
program that is well adapted to the use in a high per- 
formance computing environment. In this section some 
remarks on the computational effects of the choice of the 
underlying coupling scheme shall be given. In general, as 
long as a proper scaling is achieved the most intuitive way 
to speed up calculations using high performance comput- 
ers is to distribute calculations among many processing 
units. Nevertheless, as will be seen below the right choice 
of initial parameters like the coupling scheme can help to 
ease the problem of calculating energy spectra and ther- 
modynamic properties of magnetic molecules. 

Since several terms appear in this section which might 
be unknown to the reader, their particular meaning as 
well as related aspects shall be discussed first. The term 
computation time refers to the cumulative time that is 
needed in order to perform a certain number of fioating 
point operations (FLOPs). The execution time refers to 
the runtime of the considered part of the program. As- 
suming a parallel execution of the program with optimal 
performance, the computation time remains unchanged 
although the operations are performed in parallel. A re- 
duction of the computation time can be achieved by re- 
ducing the number of FLOPs that have to be performed. 
In contrast to this, the execution time usually decreases 
with increasing number of processing units. The change 
of the execution time as a function of the used processing 
units is called scaling behavior. It is referred to as opti- 
mal if the execution time is divided by two whenever the 
number of used processing units is doubled. The speed 



up Sp using p processing units is defined as 

n ^1 



(A34) 



where Ti and Tp refer to the execution times of the se- 
quential and the parallelized algorithm, respectively. An 
optimal speed up corresponds to Sp = p. The optimal 
speed up can be achieved if the whole source code can 
be parallelized without dependencies between the opera- 
tions which are executed in parallel. Practically, sequen- 
tial regions and communication between the processing 
units often limit the speed up to a value that is lower 
than the number of used processing units. 

In general, the computational realization of the pre- 
sented framework can be divided into two completely in- 
dependent parts. On the one hand a matrix represen- 
tation of the Hamiltonian is set up with the help of the 
irreducible tensor operator approach. On the other hand 
this matrix or independent blocks of it are diagonalized, 
i.e. the eigenvalues and eigenvectors are determined nu- 
merically. 

If point-group symmetries are used, a very decisive role 
concerning the computation time is played by the con- 
struction of symmetrized basis states. These functions 
appear as linear combinations of the initial basis states. 
The weight of the states that are included in these linear 
combinations is determined by general recoupling coeffi- 
cients. If a coupling scheme can be found that is invariant 
under all point-group operations, the number of states 
that contribute to a linear combination representing a 
symmetrized basis state is minimized. As already men- 
tioned, a reduction of computation time is achieved by 
choosing a coupling scheme that minimizes the number 
of appearing summation indices and Wigner-6J symbols. 

The construction of symmetrized basis states plays an 
important role for extending the limits of numerical exact 
diagonalization with the help of the concepts presented 
in this work. Whenever the dimensions of the appearing 
matrices are to be reduced by the incorporation of point- 
group symmetries, a certain amount of additional com- 
putation time has to be spent on the construction of sym- 
metrized basis states. Since the construction procedure 
cannot easily be parallelized, the use of more processing 
units within this particular region does not always lead 
to the desired reduction of execution time. In any case, 
one has to ensure that the recoupling formulas, which 
determine the number of performed FLOPs, are the sim- 
plest in order to reduce computation time. As already 
mentioned, this can be achieved by choosing a coupling 
scheme that is invariant under the operations of the as- 
sumed point-group. The resulting recoupling formulas 
do then not contain Wigner symbols and are optimal. 

In Fig. [19] the execution times for the determination of 
the energy eigenvalues of a cuboctahedron with 5 = 3/2 
in the subspace = 0, M = 0, Ai) are shown in depen- 
dence on the chosen coupling scheme. The used point- 
group symmetry has been D2 and execution times are 
given for the the choice of a completely invariant and a 
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FIG. 19: Scaling of OpenMP parallelized regions. Exemplary 
calculations have been performed using different coupling 
schemes (cuboctahedron s = 3/2, D2, l-i{S = 0, M = 0, Ai)). 
The crosses refer to an invariant coupling scheme while the 
circles refer to a non-invariant coupling scheme. The black 
arrow indicates the different execution times for the set up of 
the matrices and the dashed lines show the optimal scaling 
behavior. 



non-invariant coupling scheme, respectively. Comparing 
the scaling behavior, one can see that the performance 
of both calculations is limited by the construction of the 
symmetrized basis states. Furthermore, it becomes obvi- 
ous that the set up of the matrices is heavily influenced 
by the particular form of the symmetrized basis states. In 
the case of the non-invariant coupling scheme the set up 
of the matrix has been much slower than in the case of the 
invariant coupling scheme because the symmetrized basis 
states involve more states of the initial (vector-coupling) 
basis. 



Appendix B: The calculation of general recoupling 
coefficients 



In this section graph-theoretical considerations are pre- 
sented that allow to determine the action of point-group 
operations on vector-coupling states. It is shown how a 
general recoupling formula can be developed from map- 
ping general Wigner coefficients on binary trees or Yutsis 
graphs. 

In general, the problem of calculating recoupling coef- 
ficients has to be divided into two parts. The first part 
is the generation of a formula that describes the transi- 
tion between two different coupling schemes. The second 
- and much easier - part is the evaluation of a given 
formula using a specific set of quantum numbers. This 
section is exclusively focused on the first part, i.e. the 
generation of a recoupling formula that links different 
coupling schemes in systems with an arbitrary number 
of participating spins which turns out to be more diffi- 
cult. 



1. Binary trees 

In the literature one can find successful implementa- 
tions that deal with the generation of formulas for gen- 
eral recoupling coefficients which only involve a series of 
phase factors and Wigner-6J symbols .'^^'Sll xhe most in- 
tuitive way of generating a recoupling formula is to oper- 
ate on so-called binary trees. ^^^ ^^^ The correspondence 
between a binary tree and a given coupling scheme is 
rather obvious. Each coupling of two spins Sa and 55 to 
a compound spin Sc forms a triad that corresponds to 
an elementary binary tree shown in Fig. 20 This tree 



is composed of only three angular momenta and can si- 
multaneously be seen as representing a Clebsch-Gordan 
coefficient. From such elementary trees a binary tree can 
be built up step-by-step that represents the chosen cou- 
pling scheme. The tree constructed this way then con- 
tains all Clebsch-Gordan coefficients that result from the 
decomposition of a vector-coupling state belonging to the 
partic ular coupling scheme into product states (cf. Eqs. 
( [A13l ) and (|A14|). 



FIG. 20: Binary tree that corresponds to a single coupling of 
two spms Sa and Sh to a compound spin Sc- 

Operating on binary trees in order to generate a for- 
mula for a general recoupling coefficient directly leads to 
the procedure that limits the resulting expressions to 6J 
symbols. Here, the generation of a formula for the recou- 
pling coefficient {siS 2Sis^ S2SaSM\s^S4^Si' S1S2' S2SM) ^ 
that appears in Eq. (A33), shall be presented. In this 



case, generating a recoupling formula corresponds to find- 
ing a transition between the binary trees that are shown 
in Fig. [21] In other words, one transforms the set of 
Clebsch-Gordan coefficients that is represented by the 
initial tree (l.h.s. of Fig. 21 ) to the set that is repre- 



sented by the targeted tree (i\h.s. of Fig. 21). 



Following the usual graph-theoretical name convention 
the single-spin quantum numbers are referred to as leave 
nodes while the intermediate spin quantum numbers are 
called coupled nodes. The total-spin quantum number 
appears as a coupled node of a special kind and is called 
root 

There are in general two types of operations that have 
to be performed in a certain manner in order to yield the 
desired form of the recoupling coefficient. These oper- 
ations are shown in Fig. [22] and are called an exchange 
operation and a flop operation. Both operations are only 
performed on subtrees of the initial tree, thus only lead- 
ing to changes in the particular subtree while leaving the 
rest of the tree unchanged. With every operation a cer- 
tain contribution to the recoupling formula is obtained. 

An exchange operation refers to a recoupling coefficient 
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FIG. 21: Transition between two binary trees that has to 
be performed in order to calculate the recoupling coefficient 

(siS2aS'iS3aS'2S4aS'M|s3S4aS'i/ SiaS'2/ S2aS'M) . 



exchange 



(a) Exchange 




FIG. 22: Operations on a binary tree that lead to a formula 
for a general recoupling coefficient. 



that appears when considering the recoupling of two spins 
Sa and S5 within a single triad. Obviously, the only way of 
recoupling these spins is to perform an exchange between 
them. The effect of this operation can easily be derived 
from unfolding the vector-coupling states \saSbSc) and 



\si)SaSc) in terms of product states according to Eq. ( A4). 
A state of the form \saSbSc) can be written as a state of 
the form Is^s^Sc) by using the symmet ry p roperty of the 
Clebsch-Gordan coefficients from Eq. (A5), leading to a 
recoupling coefficient 



(Bl) 



In analogy, a flop operation refers to the recoupling 
of three spins 5^, 55, and Sc- Denoting the intermediate 
spin by Sd and the total spin by 5/, a successive cou- 
pling scheme would lead to states that can be written as 
\saSbSdScSf). However, a second coupling scheme can be 
designed that results in states of the form \saSbScSeSf) . 
By definition a transition between these coupling schemes 
is described by a Wigner-6J symbol resulting in 



{SaSbSdScSf\SaSbScSeSf) 



Sa Sb Sd 
: Sf ^ 



(B2) 



It has to be mentioned that the flop operation shown in 
Fig. |22(b)| was assumed to create a node, i.e. a spin 
quantum number that already exists in the targeted cou- 
pling scheme, namely Sg. Whenever a flop operation is 
performed that creates a node which is unknown in the 
targeted coupling scheme, a summation variable has to 
be introduced within the resulting contribution to the re- 
coupling formula. This summation variable is completely 
determined by the symmetry of the appearing Wigner-6 J 
coefficient. The contribution resulting from a flop opera- 
tion that creates an unknown node s^' within the binary 
tree would look like 



(-1)^^^-+^^+^^+^/) V V(2Srf + l)(2Se/ +1) 



Sa Sb Sd 



Sq S Sq 



T he des ired formula for the recoupling coefficient from 
Eq. (A33) is now obtained by performing a proper se- 



quence of exchange and flop operations. This sequence is 
in detail displayed in Fig. |23| It is not the only possible 
sequence of operations on the binary tree that leads to 
a recoupling formula for the discussed transition. How- 
ever, in this simple case the displayed sequence leads to 
an optimal formula minimizing the number of resulting 
Wigner-6J symbols. 




(e) Exchange 

FIG. 23: Sequence of operations leading to the recoupling 
coefficient (siS2aS'iS3aS'2S4aS'M|s3S4aS'i/Si5'2/S2aS'M). 
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The result of the operations then reads 

{SiS2SiS3S2S4SM\s3SiSvSiS2'S2SM) = 

{23 r + 1)(2:S2 + l)i2Si + 1)(2:S2' + 1) 



x(-l) 



51+S2+S3+S4+35 ) Si Ss ^2 I ) S2 Si Si 
S4 S Si' I I Si' S S2' 



(B3) 



This simple form was only re ach ed beca use the flop oper- 
ations shown in Figs. 23(a)| and 23(c) create nodes that 
already exist in the final coupling scheme, i.e. Sy and 

S2'. 

As long as such simple recoupling coefficients are con- 
sidered, the process of determining a proper sequence 
leading from the initial coupling scheme to the targeted 
one can be easily done by hand and does not need any 
automatization. Nevertheless, more sophisticated prob- 
lems result in transitions between binary trees that can- 
not easily be treated. Then, it becomes necessary to 
set up an algorithm that automatically creates a proper 
- ideally optimal - sec |uence. Regarding binary trees, 
known implementations'^^ of such algorithms can be 
seen as trial- and-error procedures. 

By performing a subsequence of operations, initially 
containing only one exchange or flop operation, one tries 
to find a tree containing a node that is known in the 
targeted coupling scheme. Whenever it is impossible to 
find such a tree with the given number of operations in 
the subsequence, the number of considered operations is 
increased, i.e. the subsequence is extended. A successful 
implementation of this procedure leads to a stepwise cre- 
ation of the targeted tree in which each step is guaranteed 
to be performed with the smallest (overall) number of ex- 
change and fiop operations. However, the minimization 
of the number of operations within the performed sub- 
sequences does not assure that the resulting recoupling 
formula is optimal. In general, a recoupling formula is 
optimal if the number of occurring summation variables 
and Wigner-6J coefficients is minimal. Since summation 
variables and 6J symbols are only introduced by ffop op- 
erations, generating an improved recoupling formula di- 
rectly corresponds to reducing the number of performed 
flops. 



of generating a recoupling formula was solved with the 
help of Yutsis graphs. This procedure, providing a tech- 
nically more diflicult, but at the same time theoretically 
more transparent way of generating an improved recou- 
pling formula, shall be reviewed in this section. 

The creation of Yutsis graphs is a straightforward task 
starting from the background given in App. |B 1[ In order 
to understand how these graphs evolve, an explanation 
of how to construct a Yutsis graph shall be given here. 
Additionally, the reduction of such a graph leading to 
an improved recoupling formula will be discussed briefly. 
The interested reader will flnd a deeper and more theo- 
retical investigation of general features of Yutsis graphs 
in the literature. ""^^^ 



As already discussed in Sec. A 1 an expression for a 
recoupling coeflicient in terms of Clebsch-Gordan coefli- 
cients or Wigner-3J symbols can be found by decompos- 
ing the bra and the ket states into sums of product states. 
In order to clarify this procedure, the recoupling coefli- 
cient {siS2SiS2,S2S/iSM\s2,s/iSi' S1S2' S2SM) shall be dis- 
cussed in detail as an example. The decomposition of the 
bra state was already done in App. |A 5| and is given 
in terms of Clebsch-Gordan coeflicients in Eq. (A31). 



Replacing the Clebsch-Gordan coeflicients by Wigner-3J 
symbols yields 



\siS2SiS^S2Sa S M) 



{mi} [{Mi} 

^2 S4 S 



Si S2 ^1 I f ^1 ^3 S2 

mi m2 -'Mil I Ml ms -M2, 



X I ^ I ^ \mi 171217131714) , 

M2 m4 -Mj 



(B4) 



where C{a) contains the square roots as well as the phase 
factors that appear when transforming Clebsch-Gordan 
coeflicients into Wigner symbols. The multiple sums run 
over all single-spin magnetic quantum numbers rrii and 
the magnetic quantum numbers Mi of the intermediate 
spins. The curly brackets are reminiscent of a generalized 
Wigner coeflicient. The same decomposition yields for 
the ket state 



2. Graph theoretical solution - Yutsis graphs 

As it was shown in the last section, operating on binary 
trees in order to generate a recoupling formula involving 
only phase factors, square roots, and Wigner-6J symbols 
already leads to a simple and successful procedure. How- 
ever, the process of determining an optimized sequence 
of operations remains concealed. In order to improve this 
process and thus improving the recoupling formula, ideas 
resulting from more advanced graph-theoretical consid- 
erations can be applied. In Refs. l41|112ffFl4i the problem 



\s3SaSi'SiS2^S2 S M) 



eIe 



S3 54 I / Sv si S2' 

ms m4 -Ml/ / I Ml/ mi -M2' 



{m^} (^{M,/} 

X I ^'^ ^ \ ) \m1m2m3m4) . 

M2/ m2 -M J 



(B5) 



From Eqs. (B4) and (B5) one immediately flnds the 
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following expression for the recoupling coefficient: 

{S1S2S1SSS2S4, S M\ssS4^Sr 8182^82 S M) 

X (mi m2 ms m4\m[ m'2 m'^ 777,4) 

{mi},{m'.} 



(B6) 



The abbreviations and stand for the generalized 
Wigner coefficients that depend on the quantum num- 
bers of the underlying coupling schemes which the sets a 
and P refer to. Th e de rived expression for the recoupling 
coefficient in Eq. (B6) is still involving magnetic quan- 
tum numbers. Since a recoupling coefficient is in general 
independent of any magnetic quantum number, a sim- 
plification can be found that only involves spin quantum 
numbers. Such a simplification was already found in the 
former section with the help of binary trees and will now 
be discussed on the basis of Yutsis graphs. 

The main idea of generating a recoupling formula with 
the help of Yutsis graphs is - in a first step - to set up 
a graphical representation of gen eraliz ed W igner coeffi- 
cients as they appear in Eqs. (B4) and (B5). Afterwards, 



these graphs are joined in order to build up a Yutsis graph 
that represents the recoupling coefficient. A simplifica- 
tion of the constructed Yutsis graph according to special 
operations then leads to the desired formula that is inde- 
pendent of magnetic quantum numbers. 



FIG. 24: Two diagrams of the same Wigner-3J symbol with 
different sign of the nodes. The sign is related to the cyclic 
ordering of the lines. 

The building blocks of Yutsis graphs are diagrammatic 



representations of Wigner-3J symbols. Figure 24 shows 
two diagrammatic representations of the same Wigner-3J 
symbol 



^rua rrib rricj 

Such a representation consists of three lines and one node. 
With every spin quantum number in the Wigner symbol 
a line is identified. The three lines are connected by 
the node. The node is labeled with a (+) or (— ) sign 
while the lines are characterized by the direction they 
are pointing in. The (— ) sign denotes a clockwise orienta- 
tion of the spin quantum numbers within the correspond- 
ing Wigner-3J symbol (cf. Fig. 24 r.h.s.) whereas the 



l.h.s.). The free ends of the lines represent the projections 
of the spin quantum numbers, i.e. the magnetic quantum 
numbers m^, 777,5, and 777,c. If a line leads away from the 
node, the corresponding magnetic quantum number ap- 
pears with a positive sign in the Wigner symbol, whereas 
it appears with a negative sign if the line is directed to- 
wards the node. 

It is obvious that any operation that changes the di- 
agrams of the Wigner symbol in Fig. 24 will lead to a 



Wigner symbol that differs from the original one. Chang- 
ing the sign of the node or simultaneously changing the 
directions of all lines results in a factor that can be ob- 
tained from the symmetry properties of the Wigner-3J 
symbols described in Sec. |A 1| The change of the sign 



corresponds to an uneven permutation of spins within 
the Wigner-3J symbol (cf. Eq. (A7)) while the change 
of all directions of the lines corresponds to multiplying 
the lower row of the original Wigner-3J symbol by — 1 
(cf. Eq. (A8)). Both operations result in a phase fac- 
tor of (— ij^^^f'+^c whereas any rotation of the diagram 
has no effect on the Wigner-3J symbol since the ordering 
remains unchanged. 




FIG. 25: Contracting the open ends of two lines labeled by 
Sa that leads to a summation over the magnetic quantum 
number rria of the joined lines. 

In order to construct a graph that represents a gen- 
eralized Wigner coefficient, another operation has to be 
introduced. As shown in Fig. [25j the diagrams of two 
Wigner symbols can be contracted if two lines exist that 
are labeled by the same quantum number and point into 
the same direction. The resulting graph then represents 
a summation over the corresponding magnetic quantum 
number given by 




-rua rud ruej 



(B7) 



Figures [26(a)] and [26(b)] show graphical representations 
of the generalized Wigner coefficients as found in Eqs. 
(B4) and (B5). They are easily constructed following the 



(-h) sign indicates an anticlockwise ordering (cf. Fig. 24 



conventions introduced above. 

The arrangement of the diagrams representing general- 
ized Wigner coefficients is chosen in such a way as to ease 
the contraction of both graphs and was proposed in Ref. 
11081 The graph representing the left hand side of the re- 
coupling coefficient contains only negative nodes with the 
spins being ordered clockwise around these nodes. In the 
graph belonging to the right hand side of the recoupling 
coefficient the spins are ordered anticlockwise around the 
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(a) Initial coupling scheme 




(b) Targeted coupling scheme 




(c) Yutsis graph 

FIG. 26: Graphical representation of the generalized Wigner 
coefficients for the coupling schemes contained in the re- 
coupling coefficient {siS2SiS3S2S4SM\s3S4Si' S1S2' S2SM) as 
well as the resulting Yutsis graph. 





(a) Reduction of a 2-cycle 




(b) Reduction of a 3-cycle 







(c) Reduction of a 4-cycle 



w — ^ 



+) 



(d) Interchange 

FIG. 27: Operations on a Yutsis graph and resulting contri- 
butions to the recoupling formula. 



nodes tfiat exclusively have a positiv e sign . The direc- 
tions of the lines are chosen in Fig. |26(a) as to match 
the conventional form of Wigner-3J symbols in general- 
ized Wigner coefficients. In Fig. |26(b)| they are chosen in 
the opposite direction in order to compensate the phase 
factors that result from the positive signs of the nodes. 

In Fig. |26(c)| both dia grams of the Wigner coeffi- 
cients are contracted. The contraction corresponds to 
a summation over the magnetic quantum numbers m^, 
i = 1, ... ,4, as well as over M. The summations over 
the magnetic quantum numbers of the intermediate spins 
{Mi} and {M^/} has already been included in the rep- 



resentations of the coupling schemes in Figs. 26(a) and 
|26(b)[ respectively. The resulting Yutsis graph then rep- 
resents - apart from square roots and phas e fa ctors - the 
recoupling coefficient that is given in Eq. ( B6 ) . 



In general, a recoupling coefficient is equal to a Yut- 
sis graph, which is constructed according to the above 
rules, times an additional factor. This factor contains 
phase factors and square roots that emerge from express- 
ing Clebsch-Gordan coefficients in terms of Wigner-3J 
symbols and additional contributions from initializing 
the graphical representation of the general Wigner co- 
efficients (see Figs. |26(a)| and [26(b)D . It can be written 
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N 



> 4 




^Sh Sf Se 




^ Sd Sc Se 



TABLE III: Contributions F to the recoupling formula result- 
ing from the reduction of a N-cyc\e (cf. Fig. 27). The case 
N > 4 refers to an interchange operation which is used in 
order to express larger cycles in terms of 2-, 3-, and 4-cycles. 



as the signs of the nodes have to be considered carefully. 
Within the Yutsis graph that is supposed to be reduced 
these directions and signs eventually have to be changed 
in order to match the constellation of the edges and the 
nodes shown in Fig. 27 In this process the change of the 



direction of an edge 5, i.e. a contracted line, contributes 
with ( — 1)^^ to the phase of the recoupling coefficient. 

The contributions to the recoupling formula arising 
from reducing the Yutsis graph according to the above 
mentioned operations are discussed in detail in Ref. 11081 
and shall not be further explained here for the sake of 
brevity. 




FIG. 28: Graphical representation of a triangular delta. 



as 



(-1) 



'N-2 



-.1/2 



n (25^ 



1) (25, 



(B8) 

Here N is the number of single spins in the system under 
consideration and the sum in the exponent is running 
over all intermediate spin quantum numbers Si' of the 
targeted coupling scheme. S represents the sum of the 
so-called first coupled angular momenta, i.e. the sum of 
those spins that appear in the bra-ket notation of the 
recoupling coefficient in the first position of each cou- 
pling triad. In the case of the recoupling coefficient 
{siS2SiSsS2S4SM\ssS4SifSiS2'S2SM), the sum of the 
first coupled angular momenta is given by 

5 = 51 + S^i + ^2 + 53 + 5^1' + S2' . 

The formula for the recoupling coefficient is now ob- 
tained by a successive reduction of cycles that appear 
in the graph. A cycle refers to a loop that connects a 
certain number of nodes. Depending on the number of 
connected nodes, different operations exist that reduce 



the graph. Figure |27| shows the operations leading to a 
reduction of 2-, 3-, and 4-cycles. Additionally, an inter- 
change operation is shown that can be used in order to 
express cycles which cannot be reduced immediately, i.e. 
cycles with more than four nodes, in terms of 2-, 3-, and 
4-cycles. The contributions F to the recoupling formula 
resulting from the shown operations are listed in Tab. |III| 
The values for the quan tum n umb ers k r elated to the re- 
labeled edges in Figs. |27(c) and |27(d)| are determined 
by the symmetry properties of the Wigner-6J symbols 
within these contributions. Whenever performing reduc- 
tions on a Yutsis graph, the directions of the edges as well 



A Yutsis graph is said to be reduced whenever a graphi- 
cal representation is obtained that corresponds to the one 
in Fig. 
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This representation is called a triangular delta 
and gives a fact or 1, if 5^, S5, and Sc satisfy the triangular 
condition (Eq. (Al)), and a factor other wise. 

Coming back to the example of calculating the recou- 
pling coefficient {siS2SiS sS2S4S M\ssS4Si'SiS2'S2SM) 
that is displayed in Fig. 26(c)[ one immediately finds 
two 3-cycles which can be reduced in order to generate a 
recoupling formula: S1-S2-S2' and 83-84-82' As a result 
of this reduction, the recoupling formula contains apart 
from phase factors and square roots two Wigner-6 J sym- 



bols as in Eq. (B3). 





FIG. 29: Reduc tion of the outer triangle of the Yutsis graph 
from Fig. 



26(c; 



However, as already mentioned in App. |B 1| there is 
in general more than one possibility of reducing a graph. 



Figure [29] shows a possible first step that reduces the 
outer triangle spanned by the edges si, 52, and 82'- For 
the sake of clarity the signs of the nodes as well as the 
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directions of the edges are omitted since they only con- 
tribute to the phase of the recouphng coefficient. As 
one can see after reducing the outer triangle, four dif- 
ferent triangles appear. The reduction of each triangle 
would then lead to a triangular delta and therefore to 
a completed recoupling formula. These formulas look 
slightly different but, can of course be transformed into 
each other. One possible expression for the recoupling 
coefficient is given by 

{siS2SiSsS2S4SM\ssS4SrSiS2'S2SM) = 

Si— S2— S3— S4 + 5'i— S'i/+2S'2/ 



(251. + 1){2S2' + l)(25i + 1)(252 + 1) (B9) 



S Si Sv \ S S2 

, si S2' 82] \Si 53 

So far, this section has dealt with the construction of a 
Yutsis graph and those operations that reduce this graph 
to a triangular delta. In principle, one could generate re- 
coupling formulas and calculate general recoupling coef- 
ficients with this information. In the discussed example 
of Fig. 26(c)| only triangles appear leading to an easy re- 
duction that contributes two Wigner-6J symbols to the 



recoupling formula. However, in larger systems with high 
symmetry often more complicated recoupling coefficients 
have to be calculated. In order to minimize the computa- 
tional effort it is desirable to generate a formula that con- 
tains as few Wigner-6J symbols and summation indices 
as possible. Again, they result from triangles, squares, 
and cycles of higher order. 

The most intuitive way of generating an improved re- 
coupling formula is to reduce the smallest cycles ffist. 
This idea was implemented in Refs. l41|112|lT3l and al- 
ready yields considerably improved formulas in compar- 
ison to the use of a trial- and-error technique. However, 
these formulas can be further improved by using a more 
sophisticated strategy of reducing cycles.*^^ 

Summarizing this section, one can say that with the 
help of graph-theoretical methods the effect of general 
point-group operations on vector-coupling states can be 
determined. Once this effect is known, the eigenstates 
of the system under consideration can be labeled with 
respect to irreducible representations of the point group. 
This characterization not only reduces the dimensions of 
the Hamilton matrices saving hardware resources, but it 
also provides deeper insight into the physics of the system 
arising from its geometry. 
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